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Hoofdstuk 1; Definities, notaties en afkortingen, 
Wij beschouwen uitsluitend reële konstanten, veranderlijken en funktiea, 


Notatie: 
De intervallen a < x<b, a Sx<%benasxs»b schrijven wij resp. als 


(a,b), (a,b) en (a,b). 


Definitie (1,1 


y= y(x), x € (a,b) betekent: y is een funktie van de onafhankelijke verander=- 


lijke x gedefinieerd op het interval a < x < be 


Notatie: 
n 
dy lm) n= 1,2,3sse 
dx 
ä a 
meestal schrijven wij LL ves Zl YI, 
dx dx 


Definitie WPVRE 


Een n-de orde differentiaalvergelijking is een vergelijking van de vorm: 


(a), mi 


Flater aast 


(n) 


waarin F een funktie is van de veranderlijken x,y,y'se«:,Y . 


Definitie (1.3): 


Een funktie y = y(x), x € (a,b), die n-maal differentieerbaar is en voldoet 

aan P(x,ysy',y'se.ery 0) = O voor x € (a,b) heet een oplossing van de n-de 

orde differentiaalvergelijking. 

Wij onderscheiden: 

a) een impliciet bepaalde oplossing dew.ze een betrekking van de vorm 
G(x‚y) = 0; 

b) een expliciet bepaalde oplossing y = y(x). 


Definitie (1,4): 

De differentiaalvergelijking Plet sanar bs = OQ heet lineair als 
Part sent”) een lineaire funktie is in hoen) 

De algemene vorm waarin een n-de orde lineaire diilerentiaalvergelijking 


geschreven kan worden luidt: 


sGr se rr, Gr ef ye) ble) 


fe) Rn 


Definitie (1,5): 


Als in bovenstaande differentiaalvergelijking g(x) = O dan heet de lineaire 


differentiaalvergelijking homogeen, 


Jij kunnen de lineariteit van een differentiaalvergelijking nog op een andere 
wijze definiëren: 
Laten u=u(x) en vev(x), xE(a,b) twee willekeurige n-maal differentieerbare 


funkties zijn. 


Definitie (1,6): 


De differentiaalvergelijking heet lineair als: 


nn = 


FO err s (laf) 


Als bovendien geldt: F(O,...,0,0x)=0 dan is de differentiaalvergelijking 


homogeen. 


Opgave: 
Bewijs dat (1.1) voldoet aan de definitievergelijking (1,2) van een lineaire 


differentiaalvergelijkinge 


Definitie (1,7): 


Een differentiaalvergelijking, die in de vorm 
n n= 5 
s{ ) ER )) (1a3À 
geschreven kan worden, heet een expliciete n-de orde differentiaalvergelijkinge 


Een expliciete n-de orde differentiaalvergelijking kan door het invoeren van: 


(n-1) 
S= Ig Y = JoseeesY = Yn 
geschreven worden als het volgende stelsel van eerste orde differentiaal= 


vergelijkingen: 
dy, 
dx ° Jo 
Ee 
dx © Y3 


dy 4 5 
dx = In 


ij 





dy, 


ak S= fOIgsIgreeerI)e 





5e 


Hoofdstuk 2: Eerste orde differentiaalvergelijkingen, 
Inleiding, 


De algemene vorm van een expliciete eerste orde differentiaalvergelijking 


luidt: 


dy _ 

dx es f(x,y). (2,1) 
Wij beschouwen eerst een bijzonder geval van (Bl miele? 

ay n : 

ML f(x) | (2.2) 


f(x) is continu voor a < x < be 


De oplossing van (2,2) luidt: 


y(x) = jä f(x')ax! + o, x, € (a,b) (2,3) 


c=y(x,) is een willekeurige konstante,. 
De geometrische betekenis van (2,2) en (2,3) is als volgt in te zien, De hoek 
p die de raaklijn aan de kromme y=y(x) in het punt (x,y) met de positieve 


x-as maakt, wordt bepaald door: 
teg = f(x). (2.4) 


Aan ieder punt van het gebied D: a < x < b, =o < y < too wordt een zekere 
richting, gedefinieerd door (2,4) ‚toegevoegd, Geometrisch stelt (2,2) in dit 
gebied een richtingsveld voor, 

De kromme y=y(x) met ylx)=e is nu een oplossingskromme of integraalkromme 
als de raaklijn in ieder punt van de kromme juist de door het richtingsveld 
gegeven richting heeft, 

Daar c nog willekeurig gekozen kan worden, zijn er oneindig veel integraal- 
krommen, die voorgesteld worden door (2,3). Zij het punt (E‚n) ED. (E‚n) is 


een punt van de integraalkromme: 


DN 
y= Ì f(x')ax' + ne 
6 


Door ieder punt van D gaat één en slechts één integraalkromme, Geometrisch is 
de existentie en eenduidigheid van de oplossing van (2,2) met de beginvoor- 
waarde x=X ylxp=e in te zien, 

Welke eisen moeten wij nu aan f(x,y) in (2,1) stellen opdat (2,1) voor kek» 
y(x )=y, één en slechts één oplossing heeft? 


Zij D een gebied in het x-y vlak, 


Definitie (2,1): 
f(x,y) is op D Lipschitz-kontinu in y als er een konstante ADO bestaat zodanig 
dat voor elke Gey) (x,7g) ED geldt: 


Gys) En £(x,90) | ee Aly, nm ya |- 


Existentie-stelling, 


Zij f(x,y) gedefinieerd op D en kontinu in x en y. Er bestaat dan voor elke 


(x_‚y.)ED een oplossing y=y(x) van WY. e(x,y) zodanig dat y(x )=y_. 
os dx o o 


Eenduidigheidsstelling, 


Als de voorwaarden van de existentiestelling gelden en bovendien f(x,y) op D 


Lipschitz-kontinu is in y dan heeft de differentiaalvergelijking één en slechts 


één oplossing. 


De bewijzen van deze stellingen worden later gegeven, 


Voorbeeld: 
Gegeven de differentiaalvergelijking SL Yy met beginvoorwaarden y(0)=0, 
Yy is niet-Lipschitz=kontinu voor y=0. Wij merken op dat door (0,0) twee 


integraalkrommen gaan namelijk: y(x)z0 en y(x)=X . 


Enige klassen van expliciete eerste orde differentiaalvergelijkingen die 
integreerbaar zijn, 
I, De differentiaalvergelijking van het type 
Ef = f(x)gly). (2.5) 
Deze vergelijking heet een differentiaalvergelijking met gescheiden verander- 
lijken. f=f(x), xE(a,b) en g-gly), yEle,d) met g(y)FO op (e‚d) zijn kontinue 
funkties op (a,b) resp. (e‚d). Veronderstel dat het punt (Fo) een punt is van 


een integraalkromme y=y(x) van (2,5) zodat 


vla) = Io 
Uit (2,5) volgt: ZE sr = f(x) of na integratie: 
Jay fe 
AE = } f(x)dx . (2,6) 


vS sh m 
Dit resultaat kunnen wij ook op een andere wijze afleiden, 


Wij voeren de funkties F(x) en G(y) in gedefinieerd volgens: 


F(x) == [ f(x)dx 
) ae (21) 

y(x) oe 
erde | A 


ne s(y) 


Na differentiatie verkrijgen wij: 


dF(x) 
dx =f(x) 
dG Xx 1 dy 


dx g(y(x)) dx * 





II, 


an 


enb 





De 
Optelling levert: 
ar(x) , dely(e)) _ 
dx dx n 
d 
& [r(a) + (4) = 0 
F(x) + G(y(x)) = konst. (2,8) 
Substitutie van (2,7) in (2,8) levert met de voorwaarde y= de betrekking 
(x) …= EO 
Í AL. f(z) ar 
ij gy) xe 


waarmee de oplossing impliciet bepaald ís, 


De differentiaalvergelijking van het type: 


dy 
dx 


waarin A, B en C konstanten, kan teruggebracht worden tot een differentiaal- 


= f(Ax + By + C) (2,9) 


vergelijking met gescheiden veranderlijken, 
f(u) is een kontinue funktie op (a,b). 
Stel 
u(x) = Ax + By(x) + C (2,10) 


Ae u(x )=Ax +By +C=u €(a,b). 


Differentiatie van (2,10) naar x levert: 


du _ dE … 
Tet Bs AaB fla, 


De oplossing wordt nu impliciet bepaald doors 


ra _ 
| de = X= Xs Enk 
u’ A+B f (u) 


als de integraal in het linkerlid bestaat. 


De differentiaalvergelijking van het type: 


SE = 15) (2.12) 


kan teruggebracht worden tot een differentiaalvergelijking met gescheiden 


veranderlijken, 
f(u), uEla,b) is een kontinue funktie op (a,b). Stel 


niej = ad (2,15) 
y 
y(x) Le Jo” u(x eu, 


Differentiatie van (2,13) naar x levert: 


En = = [f(u) - ul. 


De oplossing wordt impliciet bepaald door: 


je u == In En (24) 
u f(u)-u o 


als de integraal in het linkerlid bestaat. 


De differentiaalvergelijking van het type: 


U £la)y + g(x) (2.15) 
dx 
f(x) en g(x), xE(a,b) kontinu op (a,b) is een lineaire inhomogene differen- 
tiaalvergelijking. 


Beschouw eerst de homogene vergelijking 


dy; 

ee f(x) yr (216) 
met als algemene oplossing: 

Jy = C exp(f f(x)dx) EET) 


wearin ec een willekeurige konstante is. 


Stel nu dat de oplossing van de inhomogene vergelijking geschreven kan worden 


in de vorm; 

y= elx) exp (/ f(x)dx). (2,18) 
Substitutie van (2,18) in (2,15) levert de volgende differentiaalvergelijking 
voor c(x): 

SE glx) exp (-/ f(x)ax). (2.19) 
Uit (2,18) en (2.19) volgt: 

y= exp (/ flx)dx) . Jelx) exp [-f f(x)axlax = Yp: (2,20) 


Yp heet een particuliere oplossing van de inhomogene differentiaalvergelijking. 
De algemene oplossing van (2,15) is nu 

Y= Yu * Ip fa) 

of 


y=e exp(f f(x)ax) + exp(f f(x)ax) / glwexpl-/ f(x)axlax . (A2) 


Zij xela,b). Bij iedere Ve is nu de konstante c in (2,22) te bepalen zodanig 
dat vx )=ye De methode, die hier gebruikt is om de algemene oplossing van 


(2,15) te konstrueren heet de methode van de variatie der konstanten. 


De differentiaalvergelijking van het type: 


SE = f(x)y + h(x) . (2,25) 


Deze vergelijking heet de differentiaalvergelijking van Bernoulli, 





Te 


Voor d=0 verkrijgen wij type IV en voor d=1 type I, of IV maar dan homogeen, 
Stel dus at1, afO, De differentiaalvergelijking is dan niet-lineair, 


f(x) en h(x), xE(a,b) zijn kontinu op (a,b). 


Stel 

u(x) = [y()] 74. (2,24) 
Hieruit volgt: ï 

ya) = [u(a)] T° (2,25) 


Differentiatie van (2,25) naar x levert m.b.v. (2,23) de volgende lineaire 


differentiaalvergelijking voor u(x): 


Ee = (f=a) f(x)u + (1-a) h(x). (2,26) 


Dit is een vergelijking van het type IV, 


De differentiaalvergelijking van het type: 
d 
n(x,y) dE + elx,y) = 0 (2,27) 
waarin h(x,‚y), g(x,y), hier) en Saar) op een enkelvoudig samenhangend 
gebied D van het (x,y)-vlak gedefinieerd en kontinu zijn, heet een op D 


exacte differentiaalvergelijking als geldt: 


fe) dh 
n ie Bx' (2,28) 


Op D is een funktie F(x,y) gedefinieerd, die tweemaal kontinu differentieer- 
baar is in x en y zodanig dat: 
or 


òy 
Als (2,29) geldt, is ook voldaan aan (2,28) terwijl de oplossing van (2,27) 


es h(x,y); ES es gier). (2.29) 


impliciet gegeven wordt door: 
F(x,‚y) = konstant, (E50) 


De relaties (2,29) blijven geldig als wij bij F een konstante optellen. 


stel als D het gehele (x,y)=-vlak is: 


0, 


F(0,0) 


Er geldt nu: 


U 
ij 


F(x,0) 


Bi m0) ar Bi g(&,o)at (2,51) 


4 h(o,n)an . (2,32) 


ij 
Û 


(S_ar(O 
F(0,y) J ik dn 


8, 


Meb.v. (2,31) kunnen wij F(x‚y) berekenen: 


ear) oro) | Ean | nomen 


waaruit volgt: 
y _ 
g(E,O)dk + Jd h(x, n)dn. (2:35) 
Ook geldt: 


pe ar(e rX 
F(x,‚y) = F(O,y) = 5 Eind dE = ) g(E,‚y)at 


waaruit m.b.v. (2,32) volgt: 


| 


F(x,y) 


F(x,y) Ä n(O,n)dn +  eleadat. 6de 34) 


al 


Opgave: 
Bewijs dat de uitdrukkingen (2,33) en (2,34) voor F(x‚y) identiek zijn. 


Opgave: 
Bewijs dat (2,5) een exacte differentiaalvergelijking is. 


Een klasse van niet-exacte differentiaalvergelijkingen kan teruggebracht 
worden tot de klasse van differentiaalvergelijkingen van het exacte type. 
Beschouw een funktie M(x‚y)#0 gedefinieerd op D en daarop kontinu differen=- 
tieerbaar,. Zij gegeven 
d. / 2 
n(x,‚y) ern + gx,y) = 0, (2,35) 
De differentiaalvergelijking 


M(x,y) h(z,y) PE + M(x,y) glx,y) = 0 (2.56) 


is nu exact als geldt: 


ò Mh __ 9 Mg zn 
8x Ayr * (2.37) 


Voor het bepalen van de algemene oplossing van (2255) moet nu eerst de funktie 
M(x,y), die ook wel integrerende factor heet als oplossing bepaald worden van 
de volgende eerste orde partiële differentiaalvergelijking: 


MM, (2 BEM - 
Ri Ent ged. (2,38) 


Intesratie van (2,38) is in het algemeen veel moeilijker dan integratie van 
ven verste orde gewone differentiaalvergelijking, Het vinden van een oplossing 
van (2,38) geschiedt voorlopig door proberen of raden, 

Het is vaak nuttig om na te gaan of aan (2,38) voldaan kan worden door te 


stellen: M=M(x). 


Je 


(2,58) reduceert zich dan tot de volgende lineaire gewone differentiaal 
vergelijkings: 


aM Îh dam - 
Ks En = 3) = 0, (2,40) 


10, 


Hoofdstuk 3: stelsels van eerste orde differentiaalvergelijkingen, 
Definitie LS ke 


2 
Onder een n-dimensionale vector ue(Uss Ugs sees) verstaan wij een geordende 





rij van n getallen, 


Wij spreken af dat voor n-dimensionale vectoren de volgende rekenregels 


gelden: 
2) de nulvector O is de vector Del 0Os en vr 0)4 
b) zij v een vector met kentallen CvsVgseres vj) 
U= Vu, e= v, 
| i= Tesa earn 


EN eN 


u= 0 u, = Os 


En 
c) au = (au,saug,..-,au,), waarin a een Konstante is; 


d) Ë + % = (u tv Uggs eee rU VJ 5 

e) u.ve= UIV4 + UgVg + een UoVoj 

f) de norm of lengte van een vector is: 
= 2 2 2% 
lul = lu,” + U + ese Uo 15. 


Definitie (5,2): 

Zij G een deelverzameling van de m-dimensionale ruimte Cm kan bijv Î4 A of 
n+1 zijn), 

Een vectorfunktie is een voorschrift volgens welke wij aan ieder element van 
GC een n-dimensionale vector toevoegen, Wij kunnen een vectorfunktie ook als 
volgt definiëren: 


Een vectorfunktie is een afbeelding van G in de B 


Wij beschouwen n eerste orde differentiaalvergelijkingen: 
dy, (x) 

dx 
dy, (x) 


Tax 5 £(yirIgreeeIjsX) \ 63 Ĳ 


5 LCY yrIpr ee Ins) 


dr) f ( x) | 
a En SarToreee Ijs 


d 
Wij definiëren nu de veetorfunkties y en F(y,x): 
y(x) = (yy) eyg (a), ee 047 (2)) 


EGsdlfs (yy otprererIsX)e Lo (HyrTgree erg) ere gepre ren) 





Definitie (SeSlt 


— me dy dy dy 
SE he DE ü hike nt 
Eee heeft als componenten: dr © rd RT kde sen 





De vector ge y(x)dx heeft als componenten: 


rb b b 
Of ner, | Re KALE 


Stelsel (3.1) kan nu als volgt (verkort) worden geschreven: 


—> 


U») (5.2) 
Stelling (5,1): 
Te [u v| z Val |+| (ongelijkheid van Cauchy-Schwarz) (5.53) 
ii. |usvl = lul + [vl (driehoeksongelijkheid) (3.4) 
Db. en 
Khel ) yGw)axl = » ly) |ax Cie] 
Bewijs? 


2 


(aah zo Ee = 2ab 


— ve 
2lu.vl° = ai 5 uv | = 2% u; v, 5 ut, = 
i sd d 


arie 


è 2 en 

ES = 

2 B (usvj) (uv) s B (ujv;) + Cajv‚) = 2|ul 
Ly J 1e 

zieruit volgt: 


u.vl = |ul|v/. 








av |° = 5 (uy+v‚)° = lul° + kik +2 Eu,v, 5 
a: 1 
2 WÊ + el zu,vjl = lalf + WS + lul s 
5 
2 2 > > 12 
zj + WvlS + 2lullvl = {lul + Iri}. 


zicruit volgt: 


laev| < lul + Ivl . 
-b el 
sal? | [yal - 7 [0e REAL: 


z’ 


a ed 
Se bi 2 y;G)y, (B)axdk s N Bi By, Go)r, (6) [axat 


dE 


B Pa DP … te 2 

_ Prolas | [Foll - [Wola 
2’ av av a 

neeert valbete 


PE je 
{roel sf oee geerde 
a 


12, 


Definitie GENE 


Zij G een samenhangende open verzameling, kortweg gezegd een gebied, in de 


(n+1) dimensionale 449 ores 19) ruimte, 
f(y,x) is op G Lipschitz-kontinu in y als er een konstante A bestaat zodanig 


dat voor elke (x,y*), (x,y**)EG geldt: 


EG*x) — Fy) | saly* = yer). (35.6) 
Ô E 
BEANS Bf, (y4o e+ sI0*) 
Ga na dat als de partiële afgeleiden gn (i,kel,2,e.e,n) bestaan 


0) Yx 


en kontinu zijn, F(y,x) Lipschitz-kontinu is in y als G konvex isSe 


Existentiestellin 2): a 
Zij gegeven stelsel (35,1) in de vorm MF(,x). 
Zij Fly‚x) kontinu op G: ly-y, 1%, |x-x, | <a in y en x en Lipschitz-kontinu op 


G in yen zij f(y‚x) begrensd op CG, d.w.z. F(y,‚x) [sn Dan bestaat er voor 





x=, | Saa tenminste één oplossing van 
WY rr AN > DE: 
a = f(x) met y(x) = Ye (5.7) 


g=min (asp). 


Bewijs? 
Wij zullen deze stelling bewijzen m.b.v. de methode der suksessieve approxima- 
tie, ook de iteratie-methode van Picard-Lindelöf genoemd, 


Wij schrijven de differentiaalvergelijking als integraalvergelijking: 


Q 
Xx 
le) 


ek à | Free. (5.9) 
sE | 


Definieer y 


Fo Ha) + | FF) Dee 6.9) 


Xx 
le} 


x) door: 


r(2 


Definieer y 


FO Fa) ef FED, Dae (3.10) 


XxX 
le) 


x) door: 


Algemeen: 


=l_n $ had ad lan _ÍÁ 

SD) vx) + [roe )(E) Bar ne2,3rAs… (Get) 
%o 

Wij zullen aantonen dat de rij vectorfunkties FD) voor nee naar de 


limietfunktie z(x) nadert; 














rim FP (a) « Zx), 


noo 
die oplossing is van de integraalvergelijking (3,8) en daarmee ook oplossing 


van (5,7) 
69) = Frl = [egale = Mx-x, |. 


Verder geldt: 


IA 


teld Ed 


[re , ) FPD EE) at 


Oo 


Xx 


mr za) 


HA 


rl SOE =(n-1) lat 


(kex_) 2 


(3,9), (5.10) — I5 en )| = À [ir )_5 las s MA en ë 
a 
Door middel van volledige inductie volgt nu: 


pal) _ zel sm E 


(Bett) — ) - ye) | = |x-x, |M < b. 


wij beperken ons tot het interval [ex |Saemin(asg). 


zeschouw de reeks: 
ID = Fa) + GD - Fa) + F0 - Fa) + 
Old a PPN ane Bat a PE ii 4 wor 


Fa) + 2 GD) FD). (3,12) 


n=0 


leze reeks wordt absoluut gemajoreerd door de konvergente reeks? 


5 He tiene Ha ad Ea a jn y(x) | + E (eo) 
Pe 2 3! ee TET ene LE EN . 


‘zt je analyse volgt nu dat (3.12) absoluut en gelijkmatig (uniform) konver- 


dert op [xx | <a. De partiële som van (3,12) is ye () 
lim FD) (2) 
noo 
rectaat en is gelijk aan z(x), de konvergentie is gelijkmatig (uniform) op 
x=x, | <a. | 
A ad 
zen en FO Fa) + [EG (Dar, 


Hieruit volgt wegens gelijkmatige konvergentie: 


EG) =H) + je Lin EGTE er 


bd no 
le) 


en daar f(y,x) kontinu is in y: 


B) He) + | Faam FED) Dee 
B) Fa) + FE Dee. 


Deze laatste integraalvergelijking is equivalent met: 
az(x) wd bed 
dx ” f(z(x),x) 
mea.we z(x) is oplossing van de differentiaalvergelijking. 


Opmerking: 


De existentiestelling kan ook bewezen worden zonder gebruik te maken van de 


Lipschitz=-kontinuïteit. 


Eenduidigheidsstellin „53 
Als er voldaan aan alle voorwaarden genoemd bij de existentiestelling dan is 


_ 


er één en slechts één oplossing van het stelsel waarvoor geldt x)-Io 


Bewijs: 
Stel É(x) is een andere oplossing: 
— NK ed 
EO) -Fa) rf FE Dar 3.15) 
%o 
_(n+1) — Xe (n) 
Ge) =D ie) + | EG, Da (3.14) 
x 
o 


(5.13) = |E) - ya) S Max, | 


6154610 = [EFD len EO Fo lee. 


xXx 
ie) 


Hieruit volgt: ard 
Ger) 


re) = €) | En MAL CTED 


IA 


lim Fies = €) | = |z(x) - €lx)l = 0 - z(x) = Cx). 


ne 


15. 


Definitie EEDE 


f(y,‚x) is regulier analytisch op G als F(y‚x) een willekeurig aantal keer 


kontinu differentieerbaar is naar x en y, (Ap amaide 


Stelling (S.Â): 
Als £(y,x) regulier analytisch is op G dan bestaat de oplossing van 


dy _ zr > Ne: 
WY FG,x) net Fx) = 7, (3.15) 


voor [xx | Sat en is willekeurig vaak kontinu differentieerbaar naar Xe 


Bewijs: 


Als wij kunnen aantonen dat f(y,x) Lipschitz-kontinu is in y dan kunnen wij 


met behulp van de methode van Picard-Lindelöóf het bewijs leveren, 


f(y,x) 5 (Eger erIgr oee er En neee Is) 


or, Bf, 
Ee is kontinu dus begrensd: | < M, 
al Ĳ dd 
| f, 
pn % _ == Ss —_y 
(F4 o9or ee eInr®) £GAoTge 0001000) | B dy, dy, s M|y, ye]. 


1 
Hieruit volgt: 


Ero otortgeeeertgs®) = £j(rhorderoee eerde) | = 


|E (oaorortgeeee Is) = Ei(rhordeieee erge) | + IE, foretzree Ins) » 


Lyre ee00In0) | 
+ EyAorZorgoee 09) een LATE d geen 0900) | + ese 


n 
ME ly, = vile 








[Etre err) TE (rfreeerrde)| S RA 
Er volgt nu voors 
bp ee — — 2 ad hed 
EG) - EGO] = UG) - GE + IG Fr) [PJP 
n 2, & n 
sSnrl 2 woeytl| s\aMm 2 |y.-y* 3,16 
6 Liet vil i=1 vs-rfl le 
Verder geldt: 
=e aL 
yesil = : ly, rt] 
i=1 
of 
lyj-vtl = Iy-y*). (5.11) 
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Uit (3,16) en (3,17) volgt nu: 
|f(y,x) - F(y*‚x)| SnVn M[y-y*). 
De Lipschitz-konstante A uit (3,6) is dus gelijk aan nYn M, 


Wij moeten bewijzen dat y(x) willekeurig vaak differentieerbaar is. Dit moet 


dan gelden voor iedere komponent van y(x). 


dy. 
1 d e . 5 
ae ECoYgoee 09) f‚ is kontinu differentieerbaar 
volgt: 
2 
d Ji - òf, . B hi À dy. 5 df, N de df, e 
ax dx je dy; dx òx j=1 dy; J 


Deze laatste term is wederom kontinu differentieerbaar. Met behulp van volle- 
dige inductie kan nu aangetoond worden dat y; (x) willekeurig vaak kontinu 


differentieerbaar is naar Xe 


Lineaire stelsels, 
Definitie MD 
Een lineair inhomogeen stelsel differentiaalvergelijkingen is een stelsel van 


de vorm: 
== ass (274 + a,o(x)yp EE Kan arn (2)Y, d b,(x) 


: (3,18) 


wen an) + a2(*)yz + ese ann) + b, (xx). 
Als b, (x)=0 (i=l,...‚n) heet het lineaire stelsel homogeen, 


In vectornotatie luidt (3,18): 


> 


Een = A(x)y + b (3.19) 
waarin ye (y4 sees) b=(b,(x), 0. …,b, (x)) en 
Alx) = la; 3} een nxn matrix, 


Beschouw de homogene differentiaalvergelijking 


head 


z = Alx)y. (3.20) 
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Stelling (5,5): 
y=Ö is een oplossing van (3,20). 
Als zi) (v=l,eserp) voor het interval (a,b) oplossingen zijn van (3,20) dan 


is ook voor willekeurige konstanten 8, de lineaire kombinatie: 


F0) = e, 70) G.en) 


oplossing van (3,20). 


Definitie KANDE 


De vectorfunkties Pd Pl heten op een interval a<x<b lineair 
afhankelijk als er konstanten c, (v=1,eve,p) bestaan met E le, |>0, zodanig 


dat: 
Pp 


B ce, EO re = 0, (3,22) 
vel 


Bestaan dergelijke konstanten niet, dan heten de vectorfunkties lineair onaf= 


hankelijk; of de vectorfunkties heten lineair onafhankelijk als uit (3,22) 


volgt: 
ed Velg eees De 

Stelle vand 

Laten y°\/(x) (vel,e.e,p) oplossingen zijn van (3,20) voor a<x<b. 

Als geldt: : 
P =( v) = 
El A (Cx) = 0 a<x <b 


dan geldt dit voor alle xE(a,b). 


Bewijs: 
y(x) = Ee e, ars is ook een oplossing van (3,20) met y(x )=0. y(x)=0 is 


ook oplossing van (3,20), Volgens de eenduidigheidsstelling moet nu gelden: 


P jp _ 

2 ce rr = 0, 
v 

vel 


REE alde 
L MG) (val,eee,ps Pèn) oplossingen zijn van (3.20) voor a<x<b. 


aten y 
De vectoren zu (x) zijn dan lineair afhankelijk. 


Bewijs: 
Kies x_Ela,b), 


De vectoren Dar Pa) met n kentallen zijn nu, omdat p>n lineair | 


afhankelijk, dew.ze: 
Ko, zl 5 
ie dE Cx) = 0, 


ve 
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Uit de vorige stelling volgt dans 


Pp _—> md 

DH & zh = 0. 
v 

ve Î 


Svelling (3,8)e 

Zij Ax) kontinu op a<x<b. Stelsel (3,20) heeft dan op het interval ab 
precies n lineair onafhankelijke oplossingen y Î Gat rds 

Iedere oplossing y(x) is te vormen uit een lineaire kombinatie van deze lincair 


onafhankelijke oplossingen; mea.we er bestaan konstante En zodanig dat: 


Fo) = 2 oF). 3.23) 
p=Ì 


Bewijs: 


Beschouw de n eenheidsvectorens: 


e= Liemers d) 
Sig En (De tras dd 
ain = de nk 


Deze zijn lineair onafhankelijk. Bepaal vervolgens n oplossingen AE van 
de differentiaalvergelijking, zodanig dat: 


Na) «8, (vats onvsn). 


Daar er voldaan is aan de voorwaarden van de existentie- en eenduidigheidsstel- 
ling bestaan deze oplossingen en zijn eenduidig. Zij zijn lineair onafhankelijk 


want zij zijn lineair onafhankelijk voor Keke 


Definitie ‚8 

n lineair onafhankelijke oplossingen van een lineair homogeen stelsel heten een 
integraal basis of een fundamenteel stelsel, 

Toets voor de lineaire onafhankelijkheid van n vectorfunkties, 

De n vectorfunkties y (x)...y /(x), kontinu differentieerbaar op (a,b), 


zijn lineair onafhankelijk als 


n 
DS & FP) (0) = 0 

pst P 

geen andere oplossingen heeft dan 60, p=lsseerns Dit betekent dat voor 


xEla,b) moet gelden: 


Det ly, P(x)| + 0, (3,24) 
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Het oplossen van een beginwaardeprobleem, 
Gegeven; 
dy . —> — _ (0) Dn 
Teamr vla) er |= Cgseevst) 
tij =(n) 
en het fundamenteel stelsel y° ’(x),sery\ ‘(x). 
Beschouw 
Fe ze FP). 
pe] 
Er moet nu gelden: 
n 
70) En 5 ä FP) 
p=Î 
of per komponent: 
n (») De 5 
n= Be Ji Go) Ket pn ole (3,25) 


p=| 
Daar voorwaarde (3,24) geldt, kunnen uit dit lineaire stelsel vergelijkingen de 


konstanten Be worden opgelost, 


Het oplossen van een lineair inhomogeen stelsel differentiaalvergelijkingen, 


Gegevens 


_ 


LL _ AXT + D(x) b(x) = (hb, (a)s eee bx) (3,26) 


dx 
yn) = (e ) 


en het fundamenteel stelsel St a te) van de homogene vector 
differentiaalvergelijking. 
Wij passen nu op dit stelsel de methode van variatie van konstanten toe, 


d‚w.Ze beschouw: 
n 
So) Be) FP). 
peel 
Substitutie in (3,26) levert: 


n =p) n de, (x) 
Bol) GE), 7 RE reg a ze FP) + Ea, 
Pp dx 
p= p=| p=1 P 
of 
8 ge 5 de la) A 
ee) De a EF «B 
Bee 
en daar FP) oplossingen zijn van het homogcne stelsel volgt: 


a de te) ke 
EL te e= Md (3.21) 


p=Î 


of per komponents: 


20. 


n dce_(x) 


ö 
p=1 dx 


Uit stelsel (5,28) kunnen, daar voorwaarde (3.24) geldt: 


Pa) = h(x) delmanbate (3.28) 


de (x) 
ET: Pd Lr ieertt) 
opgelost worden, 


Wij vinden dan relaties van de vorm: 


de (x) 
in = 9) p=l,eee 40; (5.29) 


waaruit on onmiddellijk bepaald kan worden, 


Definitie (5,9): 
De nxn matrix b(x), met als kolommen n lineair onafhankelijke oplossingen 


y(x), velsen van (3.20) heet de fundamentele matrix. 


Opmerking, 
Als wij de fundamentele matrix zodanig konstrueren dat geldt: bx, )=1 


(eenheidsmatrix) dan kan de oplossing van (3,20) met beginvoorwaarde 
ya) 0) geschreven worden als: 

y(x) = o(x)7t°), (3.50) 
Door substitutie van (3,30) in (3,20) volgt dat b(x) oplossing is van de 
matrixvergelijking 
d® . 
desde MeT (3.51) 


Het oplossen van een lineair inhomogeen stelsel m.b,v, de fundamentele matrix, 


Gegeven de matrix Ó(x) gedefinieerd door (3.31). Zoek naar analogie van (5.30) 


oplossingen voor (3,26) van de vorm: 
y(x) = 9x) u(x). (3.32) 


Substitutie van (3,32) in (3.26) levert: 


du, dè — > > 
btn tsÂdu+ biel 


meb.ve. (3,351) volet: 


Wel) Ez) = ©). G35) 


Integratie van (3,53) levert: 


Eee Parte Been. 


xv 
le) 


u = 


ei, 


(3,26) heeft dus als algemene oplossing: 


y(x) = (rt) ii Ea ò(x)0 | (s)D(s)äs | (5.34) 


mits de matrix 8” Î(s) bestaat (Ö(x) moet dan niet-singulier zijn). 


Definitie „10 


Gegeven: 


rd > hind 
SL = AG) Alz) = (2,70) }s y= (yere) 
x=X, heet een regulier punt als alle a, in de omgeving van x=X, in een 


Taylorreeks te ontwikkelen zijn met konvergentiestraal p>0, 


Stelling { 5 9): 


Segeven: 
dy, N 
ai È 2,397; iele 
j=1 


vC) = Yze 
In de omgeving van een regulier punt konvergeren de Taylorreeksen 


n= zr ad a)". (5.35) 


n=0 
De konvergentiestraal is niet kleiner dan de konvergentiestraal van de 


Taylorreeksen voor EL 


j 


ewijs: 
Stel xd. Ga na dat dit niet ten koste gaat van de algemeenheid van de stelling. 


| 


De Taylor-coëöfficienten voor 











me 
slale Rw x (5e 36) 
d n 
n=0 
zijns n 
De 
Gú Î 
ie Bt ax [x=o 
nu geldt: 
dy; N 
dx |x=o © ‚a 2,10), 
del 
as N d da 
De ts en (3e) 
= B $a,.(0) + Y 
2 : ij dx dx J 
dx |Xx=0 j=l = X=0 





* 
etc. 





22, 


Afschatten van (3,37) levert: 


N 


dy; 
Re la, (0) | vs |. 


dx 


ÎiA 








X=0 





Beschouw nu het majorerende stelsel: 


dn, N 

ie 8 b‚;Gn, delende {Ja d8) 
j=1 

n‚(0) = |v‚| 


met als oplossingen: 


oo 
























































no) zg, Her, (3,59) 
n=0 
n 
Re! 
ï 
n DL a Iseo 
dn, N 
Tax [xeo © „5 Pa) Is | 
J=1 
2 
d N; N dn. d b,. 
AE Or al AK 
dx X=0 d=1 d =O =O 
etc, 
Stel nu: 
| db; laPa,. 
b,.(0) > |a,.(O) |; > ' (3,40) 
1J 1 ax? |x=o | dx” |xeo 
Dan volgt: 
dn, N N i 
Di lee” Db 2 5 a Og 2 TE lean 
j=1 dl 
of 
Sn K SE rz N 
dx |x=o dx jx=0 eens 
en met volledige inductie: 
n n 
.d J; d IN 
DE < n neÎs es 
dx X=0 dx X=0 














Er volgt dus dat (3,39) machtreeks (3,36) absoluut majoreert. 


De vraag is nu of een stelsel als (3,38) bestaat. 


ì 


\ 


ee AN 
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Zij 


8 


a (2) x° konvergeert voor |xl <p. (3.41) 


2, = ij 


1 4 


n 
Dit betekent dat moet gelden (wortelkenmerk) 


le) 


le AE den (M = konstant). 
ij n 
p 
Stel nus zet 
bode Leet 1 woar Jel £p (3.42) 
ae 1-À n=0 p 
o . 
en duss 
db, . xx 
= n en 
ax X=O0 or 


Uit (3,41) volgt: 


afa..(x) 
id = nl! a, (3,43) 
n ij 
dx X=O 


Uit (3.42) en (35,43) volgt dat voldaan is aan (3,40); meaew. wij hebben het 


majorerende stelsel gekonstrueerd: 


dn, M N 
Vdx Tx E nj i=lgeee,lN (5.44) 
1m j=Î 
p 
Stel N 
5 NN. 5 Y, 
ged % 
Dan volgt: 
af _ _MN = Xy=PMN 
EE Y of Y-Y(t- 5) ' (3.45) 
p 
Substitutie van (3,45) in (3.44) levert: 
dn, 
> ® M y= M (1 HHN 
Xx Xx o X p 
1 je 
p p 
of na integraties 
Yy y 
ee _ X-PMN a 
Be 7 N oe 5 + Lv; | N id 


Hieruit volgt dat voor [x | <p elke n, in een konvergente machtreeks te ontwik= 
kelen is, Elke machtreeks voor IN majoreert absoluut de corresponderende 
Taylor-reeks voor Vie De machtreeksen voor Ys konvergeren dus absoluut terwijl 
de konvergentiestraal niet kleiner is dan pe. p is de konvergentiestraal 


waarbinnen de machtreeksontwikkelingen voor Bis konvergeren. 
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Lineaire stelsels met konstante coëfficienten, 


Wij beschouwen het volgende lineaire homogene stelsel met konstante coëfficienten: 


dy, 
ES Saa * dygtp Feen Sraa 
dy, 
dx 7 Pn1/1 * Aned2 + eee Ann’ 


of in vectorvorm: 
Ed 


dy De _> 
dx si Ay (3,46) 


Eed > 
waarin Â een matrix is met elementen as en y de vector Y=(Y41Yart 000) 
Wij zoeken naar oplossingen van de vorm: 


en | (3.41) 


ed 


y=ce 
waarin cmlessCgrers0) een konstante vector is, 
Substitutie van (3,47) in (3,46) levert: 

ln je et} | (2.48) 


waarin I de eenheidsmatrix is, 


Uitschrijven van (A-AI)c=0 levert: 


ee KEN Heee Ane, = 0 ) 

ant, + (ann=Â)C) + eee An, CG, = 0 

211 2e 2 2n n 

. (3,49) 
214 + neo de ade (annen = 0, 

(3,49) is een stelsel lineaire homogene vergelijkingen met onbekenden CqseeesCne 


Dit stelsel heeft alleen dan oplossingen die ongelijk zijn aan de nul oplossing 


als de coöfficienten-determinant gelijk is aan nul: 


(as4-)) ao ajz eer Bin 

an (a5=)) 23 «1* Bon 

ee je 0 (5,50) 
dad …esesseeneeseseee (a) 


(3,50) is een n-de graads vergelijking in À en heet de karakteristieke vergen 


lijking, 
De vergelijking heeft n wortels An p=l,eees,n, die de eigenwaarden heten. 


Wij onderscheiden twee gevallen: 
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Ï, Alle n wortels van (3,50) zijn verschillend. Voor iedere À kunnen wij stelsel 
(3.49) oplossen. De rang van de matrix (A-AI) is gelijk aan (n-1) dus voor 


hek zijn oP me ® op een multiplikatieve konstante na bepaald, 


stelling SEI 


À_x 
De vectorfunkties pl), B oe p=l,eev‚n, die oplossingen zijn van (3,46) 


vormen een fundamenteel stelsel. 


Bewijs: 


Wij moeten dus bewijzen dat de n vektorfunkties 7») lincair onafhankelijk 


zijn, Me&e we 


Eg EP) oP 5 
Uit 5 B Cc ed = O moet volgen B = O P=Î,esevne 
„i Ff P 
p= 
n À_x 
5 B <p) eP 0. 
p=1 P 
Uit herhaalde differentiatie volgt: 
n „ her ! (3,51) 
ERP 0 \ 
p= P ie é 
n À_x 
2 B, <p) kee eP = _ 
p=| J 
De matrix van Vandermonde S is gedefinieerd door: 
/] 4 ee AN 
/ A, Ap bek \ 
/ 2 2 2 1 
sl M ir ae | 
| . ik 
n=j N n-1 A n=1 4 
9 ee n 
Stelsel (3,51) is nu te schrijven als: 
jp Xx nd 
oei B, e L 0) 
Ax 
0 
Bt Be n (3,52) 
. A x . 
cn) B e £ ) 0 
n 
Daar [s |+0 volgt: 
elo) Ak 
c\P/ Bg eP „jd P=l, eeen 


Pp 


Il, 


Zó. 


en omdat z)a5 volgt Pe Delen este 
QeEs de 


Als alle eigenwaarden verschillend zijn, luidt de algemene oplossing van 

OE in 
Y= 5 a c e (353) 
pet P 

waarin on willekeurige konstanten zijn. 


Niet alle wortels van (3,50) zijn verschillend, 
Met behulp van bovenstaande methode kunnen wij slechts zoveel lineair onafhan= 


kelijke oplossingen konstrueren als er verschillende eigenwaarden zijn. 


Definihie (Aetljs 
De multipliciteit van een eigenwaarde is een getal dat aangeeft hoeveel maal 


die eigenwaarde oplossing is van de karakteristieke vergelijking. 


Stelling (3,11): 
Zij ho een eigenwaarde met multipliciteit s. 


Dan bestaan er s lineair onafhankelijke oplossingen van stelsel (35.46) van de 


gedaante: 
À x 
r) ) = Fb) (») e © h=i,eee4S (3,54) 


waarin de komponenten PP) van Bh) ) veeltermen zijn waarvan de graad 


hoogstens h=1 iS. 


Bewijs: 
Zie Kamke blz. 148, 


Stelling KPRRE 


Laten Agserer A, verschillende eigenwaarden zijn van matrix A met bijbehorende 
multipliciteiten SqreresS,e 


Laten ie 


5) . bertsacard, (3.55) 


de bij A, (p=l‚es-,r) behorende lineair onafhankelijke oplossingsvectoren 
zijn van (3,46). 
Voor P=Îysverr verkrijgen wij dus S{tSpes ets =n oplossingsvectoren (3,55). 


Deze vectoren vormen een fundamenteel stelsel. 


Bewijs: 
Zie Kamke blz, 151, 


27. 


Voorbeeld: 
Bi Be a 1 
dx ET Ay; den 4792) A= (4 _3) 


Karakteristieke vergelijking: (A41) 21 > Ag gele 
’ 


Stel 1° oplossing: y = ae” an (ava) 
(3.49) -a= a(1,2). 
Stel 2° oplossing: y = (brex)e * be (b,+b,) ë = (essen), 


Substitutie in differentiaalvergelijking levert: 
e - D= Ab - (c+Ao)x e= 0, 


Deze relatie moet gelden voor iedere waarde van x, Dit betekent: 


_— 


= Ö > (I+A)c = Ö = ce = (1,2); 


4 c 
I+A)b — b = (B, 2B-y). 


PJ 


) ce + A 
) e= { 
Als algemene oplossing vinden wij: 

y= alt,2je* + (B‚2B-ve"* + v(1,2)xe"T 
of als a+B= Y: 


y= Y(1,2e + (xr, 2x- ide”, 


De fundamentele matrix van stelsel 46). 


Zij gegevens 


_ 


d _p _—> had 
Ee = Ay; (0) =c. (5.56) 
Wij onderzoeken de vorm van de matrix b(x) in y=dc en gedefinieerd door: 
dt _ : n 8 
an AEG #0) = I. (5.51) 
De reeks: En 
Äx Ax 
It Ar tees Tar tere 
! m! 


konvergeert gelijkmatig in x op ieder interval en definieert een matrix Zx). 
De afgeleide van Z(x) naar x iss 


äz Amm-1 ii oml m1 7 


_ aal — A? 
Ex A + ee. “Cme1)! dese he dese EDE + KE AZ. 


Daar Z(O)=I volgt dat Z(x)zölx). 


De matrix &(x) schrijven wij als: 


Dx) = se (3.58) 


Het oplossen van een inhomogcen probleem m.b.v, de fundamentele matrix. 


Beschouw de beide matrices: 


d(x) Hs) en d(x+5). 


Beide voldoen aan: 


dZ 
dx 


Volgens de eenduidigheidsstelling geldt dus: 


= AZ Z(0) = d(s). 


dx) . Ps) = H(xrs) 


hetgeen overeenstemt met: 
Ax As A(x+s) 
e = @ . 


e . 
Uit (3,59) volgt: (stel x==s) 


I= Ò(x) . d(-x) 
of 
El) = Hat). 


Beschouw nu het inhomogene beginwaarde-probleem: 


War +E) Fo) - 7, 


Hiervan luidt de oplossing m.b.v. (3,34), (3.60) en (3.59): 
y= otei7 te) + hi D(x-s)b(s)äx 


Voorbeeld: 
dy, 

dx 

dy 

dx 


Dit stelsel kan geschreven worden als: 


= Yo | | y(x) « (94472) 


u 


en "4 + £(E) fi y (0) Ca,sa). 


_ 


dy en ed Kaal E (6) il 
ms AT + F(x) met A = Ae o) 


Het fundamenteel stelsel vectoren van het homogene stelsel ist 


en F(x) = (O,£(x)). 


nh = (cosx, =sinx) en FC) ) = (sinx, cosx) 


(kontroleer voorwaarde (3,24)). 


Beschouw nu: 
Fa) ze) Fa + op) Fm) 


of in komponenten: 


Yi * e, (x)eosx # e,(ax)sinx 
Yo ==c, (x)sinx + ej (x)eosx. 


Substitutie in levert: 


28, 


(3,59) 


(3,60) 


(3.61) 


(3,62) 


29, 


de de 
AE COSX + Sd sinx = O0 
de de 

Tx SRE + Tp COSX = f(x) 


waaruit volgt: 


de, x 
ene =f(x)sinx — e= | f(E)sink dk + «, 
de, x 

aen f(x)cosx — e= J f(E)cosk dl + ae 


Wij vinden dus als oplossing: 


Yy = ACOSX + apsinx + sax f(E)cosk dk — soer f (Heine at 


=A,SiNX + A,COSX + coor | f(E)cosk dk + san f s(5)esnt dé 


Jo 4 2 


of 
Jy = AjCOsx + apsinx + RN f(E)sin(x-k)aE 
(3,63) 


Yg= =A,SinX + A,COSX + Je f(E)cos(x-E)dE, 
De algemene oplossing van de homogene differentiaalvergelijking luidt kennelijk 


a 
_ cosx sinx d — 
Ee lnr coen) tg) S Hlalnlo), 


Á(x) is een matrix waarvoor geldt: 
1 OQ 
Ho) =d %. 


De kolomvectoren zijn de lineair onafhankelijke oplossingen van het homogene 
stelsel. Bij willekeurige beginvoorwaarde a, en d, is de bijbehorende oplos- 


sing onmiddellijk te berekenen, 


Toon aan dat (35,63) nu te schrijven is als: 


FEN = BN Hod + JE Blx-t) F(Har. 


50. 


Hoofdstuk 4: Differentiaalvergelijkingen van de n-de orde, 


Wij beschouwen alleen explicicte differentiaalvergelijkingen van de n-de orde: 


an d ae ar! 
= lers 9 es ) (4.1) 
n dx 2 n-1 
dx dx dx 
met beginvoorwaarde: 
ay 
XxX = Ks y(x) 5 " axt = Tnt 1” mel sag nel 


waarin Dye ee N44 konstanten zijn. 


(4.1) luidt in vectorvorm: 


hd 


dy _ RT 
ix = F(y,x) (4.2) 
waarin: 
gei rYoreeerY) F(y,‚x) == WosIzreee Inr Ir e 00990) 
ay 
J4 - Y5 Yn = ut melee sn=Î 


met beginvoorwaarde: Tx) == Onyseees nde 


De theorie voor de n-de orde differentiaalvergelijkingen volgt dus uit de 
theorie voor stelsels van n eerste orde differentiaalvergelijkingen, 
Voor toepassingen is het vaak niet erg efficiënt de n-de orde differentiaal- 


vergelijkingen om te zetten in een stelsel. 


Opgave: 


Formuleer zelf de existentie- en eenduidigheidsstelling voor (4.1). 
Formuleer ook de existentiestelling voor analytische funkties op blz. 15 


voor (4,1). 


Lineaire differentiaalvergelijkingen van de n-de orde, 


De lineaire differentiaalvergelijking van de n-de orde 


ar a! 
En: ek fs) Bal tess ÊY= g(x) (an 5) 
x dx 


is na invoering van de differentiaal operator 


n= 


1 
a + Ets) a + eee f 
‚dx dx 


L 


EN 


o 
te schrijven als: 
Ly = g(x). (4.4) 


Definitie (4,1): 


Als g(x) = O dan heet (4.3) homogeen. 


Als wij (4.3) als stelsel beschouwen, kunnen wij de stelling (3,8) op blz. 18 
direkt toepassen? 

Stellin PN 

Laten LG renesf (x) kontinue funkties zijn in x op (a,b). 

Ly=O heeft dan op het interval a<x<b precies n lineair onafhankelijke oplos- 


singen o,G) p=ls «ses Ne 


Bewijs: 
Ga zelf na, 


Wij zullen een voorwaarde afleiden waaraan de n oplossingen 9, moeten vol- 


doen. 
Stel: Î 
— adi, 
y= (Y4rIorrerrIn)  Y= Y4t et 5 Ts riaatie 
Voor het stelsel zijn: 4 
do 42” 'o 
(1) 1 1 
ba = (0, (x), Ax weep 
dx 
n=Î 
dop d Pp 
oe 02 2 2 
Fe (oale mar) (4.5) 
e dx 
. n=Î 
nh ee 
J a hk 7 
dx 


de n lineair onafhankelijke oplossingen. 
Uit 

a —(i) ze 

=Î 

volgt B‚=0 i=1 ee ‚Ne 


(4.6) kan opgevat worden als een stelsel van n lineaire homogene vergelijkingen 


1 


in B,- Âls wij voor B, uitsluitend nul oplossing willen hebben, geldt: 


Ps Pos ee PD 
do do do 
d 2 
W( oss Ppree es?) = Tax Tax . ee +O (4.7) 
n=-j n=í n=-1 
d 9 d P d @ 
ie $ was ee. ee 
axt Î at | ax? 1 


Definitie 2 
W(Pgs ess 0,) heet de determinant van Wronski of de Wronskiaan, 
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Opmerking: 
Een direkt gevolg van stelling (3,6) is dat als W=0 voor x=X E(a,b) dan is 


Wz0O op (a,b), 


Stelling (He)e 


Zij px) een oplossing van Ly=g(x). °, (x) i=Îlsseern zijn de lineair onafhan= 


kelijke oplossingen van Ly=Û. 


De algemene oplossing van Ly=g(x) luidt: 


y=) + B 2,0) (4.8) 
i=l 


waarin av, i=l,...,n konstanten zijne 


Bewijs: 


Ga zelf na, 


Het oplossen van een beginwaardeprobleem, 


Gegeven de beginvoorwaarde: 
dels Hels 


Uit (4,8) kunnen wij een lineair inhomogeen stelsel algebraische vergelijkingen 
afleiden, De algemene oplossingen hiervan Casse) bestaat dan en slechts 
dan als de coëfficienten determinant van het homogene stelsel #0, hetgeen 


overeenkomt met voorwaarde (4,7). 


Verlaging van de orde van een lineaire differentiaalvergelijking,. 

Wij beschouwen vergelijking (4,4). 

Als van de homogene differentiaalvergelijking bijv. door raden, een oplossing 
bekend is, kan de differentiaalvergelijking voor de overige (n-1) lineair 


onafhankelijke oplossingen tot de (n-1)-orde gereduceerd worden. 
Laat y=y4 (2) oplossing zijn van Ly=0, 
Beschouw voor de algemene oplossing van (4,4) de transformaties 
Ge) = 1,(x). (2). (4.9) 
n maal differentiëren van (4,9) levert: 
ger 
y'= y4"P + Oye! + ye 


y'e vet yo! + (Dry or ES ye 


sy, ©), + Ey, PD pp + wan Ey olrd) + zo), 
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Substitutie hiervan in (4,4) levert: 


olLy,) + © sn dine (ll) + zak + pl) + ey mala). 


Daar Ly,=0 volgt met p'=Y: 
pnt) p pln-?) rr sh, ' (4,10) 
n-2 o Jij \x 
Voor Y is dit een vergelijking van de orde (n-1). 


cen tweede-orde Jineaire differentiaalvergelijking,. 





De algemene vorm van een tweede-orde lineaire differentiaalvergelijking luidt: 


y" + f‚lx)y!' + £_lx)y = hÚx). (4.11) 
Stel dat Jy PN Yo twee lineair onafhankelijke oplossingen zijn van de homogene 
differentiaalvergelijking (4.11). Er geldt dan: 


ie) 


rt sE rt et Lal, vp 


0, 


ij 


sn t 
lr) 4 Eil" + EA LF 
Aftrekken der beide vergelijkingen levert: 


[ror4” Atol + £yla) (74'ro - Io'rjl = 0 


of 
en + £,Gao)w = 0 
waarin 
We Fils’ — JoJ4'e (4,12) 
De determinant van Wronski is dus gelijk aan: 
we corp(-f £, ax) 
waarin c een konstante is. 
Tevens volgt uit (4,12): 
w= EU : Ent Et . (4.13) 
De algemene vorm (4.11) is m.b.v. de Liouville-transformatie te herleiden tot: 
y' + F(x)y = g(x). (4.14) 
Stel: 
y= Pe 


Substitutie in (4,11) levert: 


Li î 71 ' = 
oP + (£ro +20)! + (op "sf PHP 'f JF = h(x). 
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Kies nu: 
PB > pen _. 
Lite dep e= exp (-À | eme. 


Voor Y verkrijgen wij nu de volgende differentiaealvergelijking: 
p'" p! 
pn En + of, + f k ik . 
Ps Í o 


Pe af Pe Xx 


Verlaging van de orde, 
Stel Y4 is een oplossing van de homogene vergelijking: 
Plaka, F(x)y, = 0, 


Beschouw de transformatie: y(x)=y, (2) o(x). 


Substitutie in (4,14) levert: 


Y4P! + 2y,'9' = glx) 


of 


A: 
ak Y4 P' = 8lx) y(x) 


waaruit volgt: 


xX 4 Xx! x 4 
P =| EN Ë y,(Helk)ar dx! + ez ee dE + Eye (4,15) 
Ot y, (6) 
Uit {A,T3) vodets 
1 1 afz] 
ne De 4 4.16 
7 w dx Li ) 
Substitutie van (4,16) in (4,15) levert: 
we ser Se N 
pe En Ex Gt) tre 
EE ie aal 
“wr y,(HelbjaË - vol Helblat + ee) Ya ia, 
of met y=I49" 
PT ] 
ge) 720), (6) - vvo (6) | g(6)aE + ey, + CI, (4.17) 


(4.17) is dus de algemene oplossing van de inhomogene differentiaalvergelijking 
(4.11) waarin 4 en Yo twee lincair onafhankelijke oplossingen zijn van de 


„ 


homogene differentiaalvergelijking, 


Konstruktie van oplossingen van een n-de orde lineaire differentiaal 
vergelijking met konstante coëfficienten, 


Wij beschouwen: 
(n) (n-1) 


y + ens, + eee ay 5 0 (4.18) 


waarin a, (1=0,.evn=1) konstanten zijn. 
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Naar analogie met stelsels stellen wijs 


À 
y= e ir (4,19) 
Er volgt: 
n n= 
N ka + ese &,= 0e (4,20) 


Deze karakteristieke vergelijking heeft precies n wortels A, (ilena 


Zijn deze wortels verschillend dan hebben wij n lineair onafhankelijke oplos= 
singen van (4,18). 
Voor het geval dat er samenvallende wortels zijn, geldt er een analoge stelling 


als bij stelsels, 


Wij vermelden eerst: 


Stelling EMED 


Als een veelterm P_(x) voor x=X, een nulpunt heeft met multipliciteit s dan 


hebben ook de afgeleiden: 

S= 
| P‚_'(x), P"(x) one p{ Dx) 
een nulpunt voor XX 


Bewijs: 
Ga zelf na, 


Stelling Lhad)t 


Zij Â cen wortel van (4,20) met multipliciteit s. Dan bestaan er s lineair 


onafhankelijke oplossingen van (4,18) van de vorm: 
Ax 


o 
Yn Pj Ge herwaar 


waarin Be) een veelterm is hoogstens van de graad h-1, 


Bewijs: a ed 
4 _ d d 
Voer in de operator: L === + a nn tes De 
ax? n= ax! o 
Er geldt nu: 
bear P_0). 


Differentiatie naar À levert: 


Àx A Ax 
gie =xe B P(A) + e P_'(2) 
ook geldt: 
4) 
Zien deteuxe 
waaruit volgt: 
Ax _ | E Tax 
xe” … L P(A) B 0 | 8 


Doorgaande differentiatie naar À levert achtereenvolgens: 
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= 
bat see (2p (A) + 2x P'(A) + P "(A) he 
n n n E 
a | 


Et 

el 

@ 
Ú 


5 2 ' n AA Ax 
= ks P(A) +BB (A) + 3x B (A) + ze (A) le 
d 
ni te Te er da) | k 
in n n 
Alle rechterleden zijn identiek nul voor Aak, daar Ao s-voudige wortel is van 
de karakteristieke vergelijking P_(A)=0. Hieruit volgt: 


Ax Ax EN Ax seÂ Ant 
Le = 0, L xe = 0, Lx e =O we Lr e = Oe 


À Xx 
Dit betekent dat derhalve e ee ook 


Ax 2 An 5-1 Ax 
xe ‚ Xx e eee X e 


oplossingen zijn van (4,20). 


Met behulp van de determinant van Wronski kan bewezen worden dat de oplossingen 


lineair onafhankelijk zijn. 


Voorbeeld: 
H t == 
y!' + ary + ay = 0 





karakteristieke vergelijking: Ne + aj) + as 0. 
1 B ik 
Ma” 2 [-2, È Ie a ta, | $ 
1 Ebel Ee > Aa, Ay.p Teëel, oplossing luidt: 
’ 
= Cc Rek + Ke 
Y ET | 95 ° 
2 2 2 é 
2) a, < Aas Stel -a, t4a=w . Oplossing luidt: 
=Ba,x-iux „ba, x+iux 
y= ec, + ej e s 
Stel ie b, + i bo; e= b, - Ì bo. volgt: 
=ba,X 
y=e Jzn, COSWX + ab, sinax} ° 


3) a, = da» Oplossing luidt: 


eN = 
-ôa,x ba,x 
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Het bepalen van oplossingen van lineaire differentiaalvergelijkingen 





in de vorm van machtreeksen, 


Met behulp van de iteratiemethode van Picard kunnen wij twee lineair onafhan- 


kelijke oplossingen van 


de 
Ba + y= 0 
dx 
konstrueren in de vorm: 
1 3 ef 1E onl 
_ n+ 
Ap EE mT ie 2n+1)! 


pn 4 pen sb Xx + ze 
Jz = 2! den 


Opgave: 


Bewijs dat deze oplossingen lineair onafhankelijk zijn. (Bepaal de determinant 


van Wronski voor een zekere waarde van Xe) 


Deze machtreeksen zijn gelijk aan de elementaire transcedente funktics 


Yy=sinx en Yo=COSX. 


Het is meestal niet mogelijk machtreeksen, die oplossingen zijn van differen- 
tiaalvergelijkingen te herleiden tot elementaire funkties, Vele van deze 
machtreeksen heef t men namen gegevens vaak worden deze machtreeksen ook ele- 


mentaire funkties genoemd, Bijvoorbeeld: 


a 
ay xy = 0 (4.21) 
_2 
dx 
met als oplossingen de funkties van Airy. 


Veronderstel: 


oo 
y= 5 a x. 


n 
=O 


Substitutie in (4,21) levert: 


a nlmtj S … a x= 0 
n n 
n=? n=0 
of ka pe 
m m 
5 A42 (m+2)(m+1)x =— E an, * =O 
m=Q m=Î 
waaruit? 


co 

m 

dat 5 he? (m+2) (m+1) — an-1l Xx = 0, 
m=Î 

Deze machtreeks is dan en slechts dan gelijk aan nul als iedere coöfficient 


gelijk nul is: 


(A8) 


IV 


a, = Oi ans (m+1) (m+2) — Aj m 


58, 


Uit deze recurrente betrekking volgt: 


a, = är Ee ag = ee amel = 0 

a, = edn a a sd a 
Bei “0 zr na 4 

en Wee en NE 

6 Bb An Bedert A 


n 1 Ì 


azn Das IN et Zn+1 Bees (3n+1) lj n=l, se. 


zijn willekeurige konstanten, 


eee 


a en a 
ie) 


1 


Als algemene oplossing vinden wij: 


oo 


jn Ee 3n+1 
o Ea xXx ta, 5 B x é (4,23) 
n=0 n=0 


y=a 


Uit (4,22) volgt: 





Ant1 _ An+3 B 1 
a ĳ__@ T_(n+2)(n+3 
n n 
B a 
el Í 
B Sa 5 {n+2){n+3) * 


Volgens het kenmerk van d'Alembert konvergeren beide machtreeksen voor alle Xa 
Konvergente machtreeksen mogen term voor term worden gedifferentieerd. Achteraf 
blijkt dat het gerechtvaardigd was (4,23) algemene oplossing te noemen van 
(4.21). Het konvergentie bewijs is in dit bijzondere geval eenvoudig te geven 


omdat (4,22) toevalligerwijze van zeer eenvoudige gedaante is, 


Algemener: 


Gegeven: 
2 
EE pla) Zraloye0 (4,24) 
dre dx 


Onder welke voorwaarden is de algemene oplossing van (4,24) te schrijven als 
een konvergente machtreeks in de omgeving van KX Wat is de konvergentiec=- 


straal? Pas stelling (3,9) toe op differentiaalvergelijking (4,24). 


Lineaire differentiaalverselijkingen van de tweede orde met singuliere punten, 





Wij beschouwen: 


2 
ay P(x) dy 

red a) gen (4,25) 
dx 


waarin P(x), Q(x) en R(x) kontinue funkties op |x|<a. 
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Definitie IE 


x=0 is voor (4,25) een singulier punt als R(0)=0 en P(O)FO Q(0)+0, 


P(x Xx 
She a = p(x) en a = q(x). 


Opmerking: 


p(x) en q(x) zijn niet kontinu voor x=0; er wordt dan niet voldaan aan de eisen 
van de existentiestelling. Wij kunnen niet verwachten oplossingen te verkrijgen 
van de differentiaalvergelijking met willekeurige beginvoorwaarden van y en 


ar voor x=0, 
dx 


Wij zullen trachten oplossingen te konstrueren op O<x{<a. Soms (naar het zal 
blijken) zullen deze geldig blijven (en dus bestaan) op OSx<a, 

Het gedrag der oplossingen voor x*0 wordt bepaald door de wijze waarop 
R(x)-0 voor x-0, 


Wij beperken ons tot een veelvuldig voorkomend type van singulier punt: 


Definitie (4,4): 
x=0 heet een gewoon singulier punt indien R(0)=0; P(O)FO Q(0)+0 terwijl 
x p(x) me q(x) konvergente Taylorreeksontwikkelingen bezitten voor [x|<o. 


Opmerking: 


Wij beschouwen een singulier punt in x=0. Dit is geen beperking, immers als 


geldt R(x)=0 voor keX, dan volgt uit de transformatie X=XtX 


R(k + Xo) = R(X) = O0 voor Xx = 0, 


Voorbeeld: 


De vergelijking van Euler van de tweede ordes 


2 . 
EE 
PE er ig E: y= 0 (4,26) 


waarin « en B konstanten zijne. 


Zoek naar oplossingen van de vorm y=X « Substitutie in (4,26) levert: 


n° + (air + B= 0. (4.27) 


Deze vergelijking heet de index-vergelijking. 

Afhankelijk van de waarde van ad en B verkrijgen wij twee reële, samenvallende 
of toegevoegd komplexe oplossingen, 

Wij zullen een ander licht werpen op de index=vergelijking. Stel in (4,26) 


x=exp(t) of In x=t. Dit betekent: 


dT  ert dy 
dx ET 
2 fsd 
dy eT2t| dy dy 
8 4 e du” 
dx dt 
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Substitutie in (4,26) levert: 


2 
— Î 
e 2t ay ae (a-1) dy + ey | = 0, (4,28) 
2 dt 
+ li 
Wij hebben (4.26) dus herleid tot een lineaire differentiaalvergelijking met 


konstante coëfficienten. De karakteristieke vergelijking luidt: 


a limilke Bs Ds (4,29) 


De oplossingen hiervan zijn gelijk aan die van de index-vergelijkinge. 
Toepassing van de theorie van de differentiaalvergelijkingen met konstante 


coBfficienten levert: 


1) (a-1)° „_ 4B > O3 Ag» Ao reëel 
Ee Agt . Ae 8 A4 . a 
dk e‚ = Cy X e5 ° 
2 
2) (a-1)" — 4B = Os A, = % 
A4 t At ie À 
1 1 Ee: 
r= 6E + Ce te = Le, + Co in ad X 
5) (a-1)° — AB < 03 A, =H + iw u‚w reëel 
A = H = iw 


2 
u 7 
Y= a ej stnot + e,cosat | „xe [eysinlwinx) + ezeos(wlnx) | é 


Definitie (4.5): 


De differentiaalvergelijking 


n dy n= arl, dy 
x rn + as X Ee tove AX HK + AIT e) (4.30) 


heet de vergelijking van Euler van de n-de orde. 


Opgave: 
Toon aan dat (4.30) meb.ve de transformatie ges te herleiden is tot een n-de 


orde differentiaalvergelijking met konstante coëfficienten. 


Bij het bepalen van oplossingen van (4,26) hebben wij ons beperkt tot het inter- 
val x>0, 

Voor x<O ontstaan moeilijkheden: wij moeten dan nagaan wat de betekenis is van 
Inx en x° als r niet geheel is, Stel nu x=-E; uit x<O volgt dan E20, Substitu- 


tie in (4,26) levert met y(x)=y(-b)=u(t): 


è 
en: + al = + Bu = O3 E > 0, 
dk 5 


Oplossingen van deze vergelijking hebben wij reeds gekonstrucerd,. Dit betekent 
dat wij in bovenstaande formules Xx door -x moeten vervangen als x<0, of om 


reSelwaardige oplossingen te verkrijgen op een willekeurig interval dat de 


41. 





oorsprong niet bevat, vervangen wij de formule op blz. 40 x door |x 


Wij resumeren: 


Stelling (A58 


Stel voor het bepalen van oplossingen van de vergelijking van Euler: 
xy + axy! + By = O 


op een willekeurig interval dat de oorsprong niet bevat: ye. 


Bereken de wortels r‚ en r, van de index-vergelijking 


Flr) e= r° + (a-t)r + Be Os 


r, Fo 
Âls rr en reëel: y=c4 |z| +c |z| . 


r 
Als r‚=r„, en reëel: y=legtejlnlx}}|x| 4 


1 2 


Als r‚ en r, complex: y=|x|*(e,sin(win|x|)+egeos(win|x[) } 


waarin Ess Fo = Ht iwe 


Oplossingen in de vorm van machtreeksen in de omgeving van een gewoon 
singulier punt, 
Stelling (4,6): 


Gegeven de differentiaalvergelijking: 
Zn 2 
Ly = xy" + x{xplx)}y' + [x°alx)}y = 0 (4.351) 


x=0 is een gewoon singulier punt, dew.ze xp(x) en Zal) zijn analytisch voor 
“=O, 


De machtreeksen: 
oo 2 co 
n n 
xplx) = BE px en xqlx) = E ax ad) 
n=0 n= 


konvergeren voor o<|x|<p. 


Verder geldt: 


P, = lim xp(x), a, = lim zal), 
XO x=0 


Laten EF, en Fo (r,=r,) de reële wortels zijn van de index-vergelijking: 
F(r) = r(r-1) + Por +a,= 0, 


(4,31) heeft dan in een der intervallen -p<x<O of O<x<p twee lincair onaf han- 


kelijke oplossingen Jy en Yo van de volgende vorm: 
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a) F4 = Fo niet geheel: 
4 oo 
yy) = |xl| k + E an (r)e"} 
n=Î 
5 a (4.35) 
2 n 
Yo lx) = Lx | ii de B a, (ro) } 
n=Î 
b) Ey = Fo? 
E15 ee n 
7) e= xl tz are} | 
n=Î 
Pé (4.34) 
Yo (x) = y(x) In |z| + |x| Î 5 blj. 
n= 
ec) FE, = Fy = N, N > 0 geheel: 
Pr oo 
y,) = | | Eg + B a{(r xr | 
L n=Í as 1 dS î 
rr en (4.35) 
2 n 
Tele) = uyrelejlaiel + kel Ae 5 eirlt be 


De coëfficienten a,(r,)s a, (ro), br), er) en de konstante a (die ook 
nul kan zijn) kunnen bepaald worden door de machtreeks=oplossingen in (4,31) 
te substitueren, 


Wij zullen afleiden dat: 
da, 
br) 7 dr 





r=r 
1 


nr) = & (eer da,] 





Eg 


a = lim (rr )ape 
Pt 


Bewijs: 
Wordt niet gegeven, Een toelichting volgt wel, 
Wij beschouwen het interval x>0 en stellen: 
lee) 
ye Ear. (4.36) 
n=0 
Substitutie van (4,36) en (4,32) in (4,31) levert: 


E+N 


RE sl a, (r+n)(r+n-1)x de wen 


ar(r-ijx + a, (r+i)rx” 


4 {np apr + px + ) Da Ex + a Catt + a_(r+n)x tE l 
o 1 n ie 1 n „j 
j= et 
+ (qa +q,x + me & kla bela on hal” l= 0 
A 44 ee En vee Lo 1 . n aaan ° 


Vermenigvuldiging van de machtreeksen en rangschikking van termen van gelijke 


machten van x levert: 
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5 r+1 
a F(r)x + la,F(r+1) + a, (pyr+ag) Je 
+ [aoF(r+2) + a (por+a,) + as(p,(r+1)+q ande nn 
2 o\Por*@p 1*P4 1 
eid + (a,F(r+n) + a, (porta) + a, (py (2-1) tan} 
+ eee A 4 {P4 (r+n=1)+a xr Hie s (4,37) 
n-1 71 1 
waarin P(r)er(r-1)+p FtAe 
In een kompakte vorm luidt den 
oo 


n- 
Ly=Lx E ax = a RIOES + 5 Fame, +5 
n=0 n=Î 


4 
za) COP le ee = 0. 
(4,38) 


De coëfficient van iedere macht van x in (4,38) moet gelijk aan nul zijn. 
Daar a+ moet gelden: 
F(r) = r(r-1) + Por +a,=0 


met als oplossingen r, en To, en verder; 
n=Î r A 
F(r+n)a, + 5 a LED ee + An-k | = 0, (4,39) 
=O 


Uit (4,39) volgt dat a, een funktie is van a ..…a,_4 Mits 


P(r+1),F(r+2),...,F(r+n) ongelijk aan nul is, 


a) Daar F(r,)=F(x,)=0 volgt dat als rj+ntr, voor nzi P(r,+n)+0. De twee 
lineair onafhankelijke oplossingen worden dus voorgesteld door (4,33) waar- 


in de konstante a, op 1 genormeerd is. 


b) FP, = Epe 


Beschouw de funktie f(r‚x)=x Le + a 5 
dat F(r+n)#0, an (r) gegeven wordt door (4,59) voor n=1l, Er geldt nu: 


n : k 
(r)x } x>0, waarin voor alle r zodanig 


oo en 5 
LP(r,x)=a P(r)x + EB Fame, +, + a CP an | ge 
Ee (4.40) 
Substitutie van (4.39) in (4.40) levert: 
L Ó(r,‚x) = a,F(r)x. (4.41) 


Nu geldt: Fjaler, 1. 
Dit betekent dat He, sx)= =y,(x) oplossing van de differentiaalvergelijking is, 


Differentiatie naar r in linker en rechterlid van (4,41) levert: 
Lndend eren) =| arl) 
Pr L o 


of s 


44, 


ôr | co zi g 

k { ee r+n „ © da, nl 
Ljx In x ja, + Da (r)x PR B mX |= 

L ( n=Î n=Î al 


Yr 2 
= 2a,(rr‚)x + a, (rr) EM Es 


Hieruit volgt dat voor r=r4? 


Es ke nia n 
y,(x)ln x + x EX 
d dr | 
n=Î r=r 


4 


oplossing is van (4,31). 
c) PR We N N > 0, geheel. 
De. uitwerking hiervan is vrij gekompliceerd. Wij zullen dit geval toelichten 


aan de hand van een voorbeeld: 


De differentiaalvergelijking van Bessel, 


De differentiaalvergelijking van Bessel luidt: 


yr + xy! + en - v)y = 0 (4,42) 
of Dn 
d | ay | 2 2, 

X Ax Le ren (x° — v°)y = 0 
waarin Vv een konstante is, 


Wij zullen voor v=0, veh en v=1 oplossingen konstrueren in het interval x>0, 


Stel 


oo 


ye E ax. (4,43) 
n=0 
Substitutie in (4,42) levert: 
co ee oo 
7] 
E a, eem) = | te E ao x= 0 
n=0 n=è 
waaruit volgt: 
2 2 
a (r = vo) = 0 
2 2 2 
a, [(er1) = vj 0 (4.44) 
2 2 E Ee 
an [{men) 4] + Bn Ô nEze2 
of - 
Ee f n pe „ii = 
L4,2 sE Ve Be Hi a 5 ns 2 
r+n)“-v 


Di an s Sg =d == eee = . 
Dit betekent: a. dg de Aon+1 0} 
Stel n=2m: 


Ì, 


Taom-2 hike: 


2m 
a ee Ae 


an leem er [(es2m) ove [Cerone ) ev] 7 
(1-)Za, 


Ë meal on Mesene len tere) | 


Ad al 


… {4Â5) 


Wij vinden dus als lineair onafhankelijke oplossingen als vO en v niet geheels 


er zo) (4,46) 
yal) ex 5 ba (4.47) 


waarin b,(r)=arn 


Verder geldt: 


B ef Am + Amv Am(m+v) 
De reeksen konvergeren voor alle x20 (a'Alembert). 


Wij merken op dat er moeilijkheden ontstaan als v een natuurlijk getal ise 


Wij onderzoeken: 
v=0s samenvallende wortels r4=Fg=0. 


(4,45) luidt: 
m m 
EE 
2m (am) (2m-2)2..2 2m) 
Een oplossing is dus: 


TT oo m 2m 
vG) =a 1 + E Ea, k a JG). (4.49) 


el 
SL mei 2 Pnt) 


(4,48) 





Deze reeks konvergeert voor x>0 en is analytisch voor x=04 

Voor het bepalen van de tweede oplossing die lineair onafhankelijk is met 
(4.49) zijn twee methoden, 

Wij kunnen de uitdrukking voor yo(x) (4.34) in de differentiaalvergelijking 
substitueren en vervolgens de coëfficienten b, bepalen. 


Wij kunnen b, ook bepalen uit: 


d - 
b‚(r) = 5e ay (zr). (4,50) 
Wij zullen deze laatste methode volgen. 
Uit (4.44) volgt: dag = 0 Uit de recurrente betrekking volgt dan ook: 
dr 


Sn d 
Wij berekenen dus alleen: Sr an (2) eet 


(4.45) luidt voor v=0t 


46, 


(17 a, 


Te m=1,2,3s eee (4.51) 
(2m+r) (2m-2+r)“.(2+r) 


dom 


Merk op dat als: 
B B B 
£(x) = (xe) (aa) Se (ama) * 
dat 


f'(x B, B, n 
f(x Xt xd xXx 











Voor (4,51) vinden wijs 


ij Ì > 
SE om”) hed eee pen EDE ds} 
a, (Fr) 2m+r 2m-2+rL ER jj 


of voor r=0: 
AEN. dE 





a! (0) = —2 


2m 


Substitutie van (4,48) in (4,52) levert: 
(-1)T a 
le) 


NE. delven (dn 54) 


‘ in 
& om) = -H 


waarin 


4 1 1 
Heh pt te (4.54) 


De tweede oplossing luidt dust 
65 (tg, et 
Jo (x) = JG) In Et Dg Hs (4,55) 
m=1 2°(m!) 


In plaats van Yolx), wordt als tweede lineair onafhankelijke oplossing een 


lineaire kombinatie van Yo lx) en JG) genomen: 

an al 

1) eG iv2E) + (ram 2)3, a) | (4.56) 
y is de konstante van Euler-Mascheroni en is gedefinieerd door: 


y= lim (HE, = In n) # 0,5772, 


noo 
Substitutie van (4,55) in (4,56) levert dan: 

mt 

me oo (-1) H 

2 2 

Leijer mD 0e 2 mx EE, (4,57) 
mei 2 (m!) 

De algemene oplossing van de vergelijking van Bessel van de orde nul voor 


x>0O luidt dans: 


y= €, JG) + rx). 


Wij merken op dat Jo (a) als x-0 terwijl rx) een logarithmische singula- 


riteit heeft voor x=0, 


47. 


II. ved. 
Zit (4,44) volgt: r4,2°2® agree eon 
7 B A1 
a (r+n)° - — '= =a nz2, (4.58) 
Af. 4; n-2 
Substitutie van r=t in (4,58) geeft: 
a = el n=ê n=è,4 6 
n ° r{n+1) anika 
of met n=2m; a 
_ Some Von= _ a (-1) do 
Aom 7 2ml2me1) © (2mr1)eml2m-1)(2m-2) 7 **° T{2m+1)! 
en dus: Ei EN: „2m = 
y(x) = XP 1 + 5 ri 
1 bee 2m+1)! | 
BE hk oo ED „em+1 + 
y,G) = X Gi = x * sinx. (4,59) 


m=0 
De Bessel-furktie van de eerste soort en van de orde & is gedefinieerd als: 
Bk 
JG) = rd sinx. 
Voor r=-k zijn de coëfficienten van a, en a, in (4,44) gelijk aan nul. a, en 


ag kunnen dus willekeurig gekozen worden. Bij ap verkrijgen wij, met behulp 
van de recurrente betrekking, AprAgser-sAon en bij as: SnBgper on donad 
In de tweede oplossing komt de logarithmische term niet voor, 


Ga na dat: En 
(-1)P a, 


Aon NE 5 n=1,2sese: 
(-1)Pa, 


dona * eni)! nh 


waaruit volgt: 


5 oo n_2n oo n_e2n+1 : 
y (x) ait E „_Ì) X Een E _Ì) XxX 7] _ a SOEX dh sinx 
2 Ee 2n)! Nan C2n+1)! RN Tx Di Ia 


=O 

(4,60) 

De Besselfunktie van de eerste soort en van de orde -& is gedefinicerd door: 
= 1ZË | 
J_40) = le COSXe 


De algemene oplossing luidt: 


y= c, 3,6) + Cp Jz). (4.61) 


III, 
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Ve, 


Uit (4.44) volgt: r, o=+t1. 


ST Sg > Sand 


De recurrente betrekking luidt voor r=1: 
—a 
n=-è 


ar hln+2) Bed hnOna we 


n=èm levert: 


+a (ih a, 


—=â 
a en et er 


Det 


e, = 
En 2m(m+1) 24 (m+1)m.m(m-1) 
Stel dd de oplossing is dan: 


ne gj 2 
y; Cx) =x & en Bs (4.62) 
m=o 2°"(m+1)!m! 


De Besselfunktie van de cerste soort en orde een is gedefinieerd door: 
3) = dy, (x). (4.65) 


(4,62) konvergeert voor alle x3 J,(x) is dus gedefinieerd voor alle xe 


Voor het bepalen van de tweede oplossing substitueren wij: 


de oo sf 
Yo (x) = a J,() In x4+x k + EE ec | > | 
nei 7 
| 
in (Adf): 
oo oo 4 
2axJ,'(x) + B (n-2)n eo, ar: e,#* | =0 met ost. (4,64) 
n=0 n=0 


Substitutie van (4,62) en (4,63) in (4,64) levert: 


ln 
+1 ie On-1 | Re ae 


és [9-eote xt B 1 (nete 
vd gef a 


= oo mn 2m+1 
| 5 (-1)"(2m+1)x E (4.65) 


== == Xx + . 
L m=Î 221) im! J 


In het rechterlid staan sleet:ts oneven machten van x3 dan volgt dus uit e‚=0; 
Edens Do ade Verder geldt: a==c ==le 
Stel in het rechterlid n=2m+í1. 


Wij verkrijgen dan de volgende betrekking: 


m 
{-1) (2m+1) me sE desa 


(eme)? | c + Cn = 
L ae ze 22 m+1) Leal 
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m=1 levert: 


è (=1)3 
Ge te, he 22 er 


ce, kan willekeurig gekozen worden. 


2 


Kies c‚ = al, Dan volgt: 
2 22 





zj s 1 1 
Cc, = 14) (HE, + Hide 
4 24 os AU J 24,21 de 


Het is eenvoudig om aan te tonen dats (dit betekent dat het vrijwel niet 
doenlijk is) 

CHI TE tE) 
a MG m=1,2, 34e. 


5 2m 
2 m!(me1)! 


2m 


De tweede oplossing luidt dus: 


=e (I) UH Hs) 7 
ene ite 2 nn OE, (4,66) 


Als tweede oplossing, de Besselfunktie van de tweede soort van de orde cen, 


wordt gewoonlijk een lineaire kombinatie van J,(x) en y(x) genomen: 
gi ä 
Y‚) ei [7720 + (y = In 2), (ze) | e 
De algemene oplossing van de vergelijking van Bessel voor vel luidt dan: 


y(x) = e, JG) + e, rx). (4.67) 


50. 


Hoofdstuk 5: Randwaardeproblemen voor tweede orde lincaire differentiaal 
vergelijkingen, 


Definitie (5,1): 
De differentiaalvergelijking 


2 
2E + ax) ze Rey = f(x), Oxi (5.1) 
dx 3 


waarin Q(x), R(x), f(x) kontinu zijn op [0,1] met randvoorwaarden 


y(o) = « yi) = pf (a‚B konstanten) (De) 
heet een randwaardeprobleem, 


Definitie (59.2): 


Het randwaardeprobleem heet homogeen als f(x)=0 op [0,1] terwijl bovendien 
geldt: y(O)=y(1)=0, 


Existentie en eenduidigheid van de oplossingen, 
De stellingen die wij hierover voor beginwaarde problemen bewezen hebben, 


gelden niet voor randwaardeproblemen, 
Voorbeeld: 
y' + ty = 0 met y(o) = 0 yi) = 1. 


Algemene oplossings: y=eCcosTx+ejsinnx 


y(o) = ce, = 0 
y(1) voor alle c‚ 


Aan de tweede randvoorwaarde kan niet worden voldaan, De oplossing bestaat 


ij 
©, 


niet. 
Wij zullen aantonen dat randwaardeprobleem (5,1), (5.2) dat dus ook het speci= 
ale geval f(x)=0 op KP impliceert, herleid kan worden tot een randwaarde- 


probleem met homogene randvoorwaarden: 


a d, 
Be + QE + R(x)y = £*(x) 


y(0) = y(1) = 0 
Beschouw een willekeurige funktie g(x)€C'[O,1 ] met &(O)=a en g(1)=B 
(bijv. g(x)=ar(B-a)x |. 
Stel: 
y(x) = nx) + glx). 
Substitutie in (5,1) levert: 


ne + Qax)! + Rlx)n= f(x) - [g" + QUx)g' + Rlx)g] = f*(x) 


51. 


met randvoorwaarden: n(O)=n(1)=0, 
Laten y;) en Jo (x) twee lineair onafhankelijke oplossingen zijn van: 


y' + Q(x)y' + R(x)y = O 653) 
y(o) = a yi) = B. 


De algemene oplossing luidt: 
j => e494 (x) ka eYolx). 
Opleggen van de randvoorwaarden levert: 


le à 


€494 (0) + e5Y,(0) ie 
e‚7 1) ai e2Yo (1) = B. 
Dit inhomogene lineaire stelsel algebraïsche vergelijkingen is eenduidig 


oplosbaar in e,‚ en e5 als: 


0 
De Zabel A + 0 (5.5) 


20 10 
Voor het homogene randwaarde probleem (a==0) heeft (5.4) met voorwaarde 


(5.5) geen andere oplossing dan de nuloplossing! 


Wij kunnen nu de volgende eenduidigheidsstelling formuleren: 


Stelling kBe de 


Gegevens 

y" + Q(x)y! + R(x)y = 0. (5.6) 
Als er een oplossing van (5,6) bestaat met y(O)=a, y(1)=B en D+O dan is er 
slechts één, 


Bewijs: 


Stel p(x) en Y(x) voldoen aan de differentiaalvergelijking en randvoorwaanden, 
Zij verders Fb 
Er volgt dans 
1 | pen __ = 
Yo" + Q)y' + Rlx)y, = O F0) = 7,1) = 0. 


Daar D#O voldoet slechts Yos waaruit volgt: g=Y, 


Samenvatting: (alternatief van Fredholm) 
DE het inhomogene probleem heeft slechts één oplossing (het homogene probleem 
heeft geen aridere dan de nuloplossing), Óf het homogene probleem heeft de 


niet-triviale oplossing. 
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Opmerking: 


Als het homogene probleem een niet-triviale oplossing bezit dan kan uit het 
voorgaande niet de existentie van de oplossing van het inhomogene probleem 
worden afgeleid (deze oplossing blijkt te kunnen bestaan onder speciale voor- 


waarden |. 
Zeker is dat als de oplossing van het inhomogene probleem bestaat, deze niet 


eenduidig is, 


Opmerking: 
Zij ox) oplossing van P'+Q(x)P+R(x)p=0 


v(O)=a op(1)-B 


en Wx) oplossing ven YV'+Q(x)V'+R(x)Vef (x) 
Y(o)=0 Y(1)=0, 


Dan voldoet y=p+Y aan: 
y'r Qla)y! + Rlx)y = f(x) 


y(0) = a (1) = B. 


Inhomogeen probleem met homogene randvoorwaarden: 
y" + Qx)y' + R(x)y = f(x) 
yo) = y(1) = 0. 


Met behulp van de Liouville transformatie kunnen wij dit probleem herleiden tot: 


(5.7) 


1" + Rlx)n= fx) 
(5.8) 
ni0) = nil) =d 
y= nexp (-& f Q(x)ax)} 
R*(x) = R(x) - À Pe) EN 


f(x) = f(x) exp (5 £ QG)adx }. 


De algemene oplossing van (5.8) luidt: (4,17) 


OE zee na Ee 
zr SS … % À [ t 
nG) =d Ent n GO) np) FENDER + Agm (a) + Apnple) (59) 
waarin 1, en 1 twce lineair onafhankelijke oplossingen zijn van de homogene 
differentiaalvergelijking (5.8). waw (ns 19) is de determinant van Wronski, 


Wij konstrueren ny, en n‚ nu zodanig dat geldt: 


n‚(0) = 0 n‚(1) +0 
n‚(0) + O 1,1) = 0, 


5de 


Opleggen van de randvoorwaarden van (5.6) levert voor (5.7): 


i 7 
no(o)l & | n,(E Dat + A, |= 0 
2 L [ 1 2 | zt 





F4 1 7 
— % = 0, 
„| Ma KEDEMCJAE + A, je 0 
Substitutie van A, en A, in (5.9) levert: 
Á 5 
1 k x 
(ed | mom Oeras + | mone. 510) 
Deze uitdrukking kunnen wij ook schrijven als: 
1 
ne) = | cx, E)en(gar (5.12) 
waarin: ) ) 
“en CE) nx E<x 
TOMA (5.13) 
d 
er nye) no (6) E>x 
G(x,E) heet de funktie van Green. 
De voorwaarden n‚ (1) + 0, n‚ (0) EO: 
Beschouw twee lineair onafhankelijke oplossingen van ip en Îles van 
n'+R (x)n=0 welke aan geen van de randvoorwaarden voldoen. 
1 en No worden nu gegeven door: 
ie n (0) „ 
GD) = HG) - Fro 1202) 
1 1 No 0 2 
5 (5.14) 
no (1) 


n‚(x) = n,lx) - CAO) nx). 

De voorwaarde n,(1)+0 en no (0) +0 betekenen voor (5.14): 
n,(1)9p(0) — 0, (0)N(1) +0. 

Deze voorwaarde komt overeen met (5,5). 


Wij zullen aan de hand van een voorbeeld een aantal eigenschappen van de 


Greense funktie afleiden. 


Beschouw: 
a 2 
oe + k°y = f(x) Oxi (15) 
dx 
y(0) = y(1) = 0 f(x) € C'[0o,1]. 


k is een konstante waarvoor geldt: 


sink + O0, 
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De oplossing luidt: 


4 
y(x) = | f(E)G(x, E)AE (5,16) 


met: : 5 
E sinkE sink(1-x) E Cx 


ksink 
Blerk) = 


Ee sinkx sink(1-É) Ebo 


ksink 


De funktie van Green wordt bepaald door de homogene differentiaalvergelijking. 


Ga na dat G(x,E) de volgende eigenschappen heeft: 


1) Gt! + k°a = 0 x tE 
2) G(x,E) is kontinu in x= E 
: _ + 
3) iin G(x,E) < tin 9G(x,E) ie AEX ) " G(x,x …) ee 
dE ôË ôk dE 

Etx Elx 

De afgeleide van G naar E is diskontinu in E=x, 
4) G(x,0) ha G(x,1) = 0 
5) G(x,E) = GLE). 


Het is mogelijk de representatie van de oplossing (5,10) met behulp van boven- 
staande eigenschappen direkt uit de differentiaalvergelijking af te leiden. 
Vermenigvuldig daartoe (5.15) links en rechts met G(x,E) en integreer van O 


tot 1: 


1 2 4 i 
hike: + Ky | G(x,E)dx = Ë f(x)G(x,E)ax. (5,17) | 
dx J 0% 
Wij beschouwen de intervallen: O0 Sx Sa <t 
ESB ExE1. 
Notatie: lim f(a) = f(E ) 
atE 
lim f(@) = f(E“), analoog voor de afgeleiden. 
BIE 
1 E kan 
hi Var? | G(x,E)dx=lim [lont LG(x,E)ax+tlim ves | G(x,E)dx, 
n L ne até dl BIE L 


Na tweemaal partiëel integreren van de eerste termen van de beide integralen 


in het rechterlid verkrijgen wij: 
. + = ús | fe der | 
J £a)ola, Bax = (6) [OEL OOK LE) |+ 9"(6) [OCE E)=CCE E) |e 


+ y(0)e'(0,E) — y(1)G'(1,E) + G(1,E)y'(1) - G(O,E)Y'(O) + 


A 
+ ER (G* + KEa)y dx + El (er + ko)y dx. 


É 


2de 


Met behulp van y(O)=y(1)=0 en de vijf eigenschappen van de Greense funktie 


verkrijgen wij: 


4 t 
y(E) = Je G(x,b)ax = E/ [7 + Ky) G(x,E)ax. (Se 18) 


Wanneer wij de laatste integraal niet opsplitsen doch direkt tweemaal partiëel 


integreren, vinden wijt 
i 2 
see [or + ora 


(Graka) is overal nul behalve in x=f, deze uitdrukking is geen funktie in 
de gebruikelijke zin. 
Er kan worden aangetoond dat rake gelijk is aan de delta-funktie van Dirac. 


Deze delta-funktie Ò(x=E) is gedefinieerd door: 
dl 
| omslaat » 9x) x € [0,1] 
o 


waarin p(x) een willekeurige kontinue funktie op (0,1) is. 


Eenvoudige eigenwaarde problemen, 


Wij beschouwen het randwaardeprobleem: 


we P(x,A)y' + Q(x, Ny = 0 
nein (5.19) 


y(0) = y(1) = 0 


waarin À een parameter is, 


Definitie (5,3): 
De waarden van À waarvoor (5,19) oplossingen heeft die ongelijk zijn aan de 
nuloplossing heten de_eigenwaarden van het probleem; de ermee corresponderende 


oplossingen heten de eigenfunkties, 


Voorbeeld 1: 


y' + Ny = 0 


y(o) = y(1) = 0, 


Algemene oplossing: 
y = Asinkx + Beos Axe. (5,20) 


Randvoorwaarden: y(O) = B = O 

y(1) = Asink = 0 

sink = 0 > A = nn n=1,2,Sge.« 
De eigenwaarden zijn X = nn. 
De eigenfunkties zijn DS Asinnnxe 


De konstanten A zijn willekeurig. 
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Stel nu A= 1, n=i,2ge ce 


Ga na dat geldt: 
4 
| Sne = 
o 


k als n=m 
O als ndn. 


Beschouw: sec [o,1] terwijl f(O)=f(1)=0. 


De eigenfunkties kunnen gebruikt worden om f(x) in een reeks te ontwikkelen: 


oo 


f(x) = ZE a). 
n=1 


Er kan bewezen worden dat deze reeks gelijkmatig konvergeert op rel. Tevens 


geldt: 


4 
a = 2 4 f(E)sinnntdt, 


n 


Voorbeeld 2: 
y' + Aye 0 y(-1) = y(1) = 0 
Randvoorwaarden opleggen aan (5,20) levert: 
AsinÀ + BeosÀ = 0 
| (Be) 


„=ÂsinÀ + Beos = O 


(5.21) is een stelsel van twee vergelijkingen met twee onbekenden A en B, deze 


verschillen van de nuloplossing als: 


sinÀ COS À 


= 0 (5.22) 





„sin cos À 





of : 
sin 2À = 0 


waaruit de eigenwaarden volgen: 
22 
nn 


A, = 4 …e 





Samenvatting: 
Uit deze voorbeelden blijkt dats 


1) er bestaan aftelbaar oneindig veel eigenwaarden; 
2) de corresponderende funkties zijn op [0,1] orthogonaal; 
3) iedere Funktie fect [o,11, f(O)=f(1)=0 kan in een reeks van eigenfunkties 


worden ontwikkeld, 


27. 


Deze eigenschappen blijken voor de volgende klasse van eigenwaarde problemen 
te gelden: 
wit + (Aplx) + Q(x))w = O 
w(0) = w(1) = 0 (5.25) 
Q(x) en plx) € C'[0,1h plx) >O op [0,1]. 
Dit eigenwaardeprobleem heet het probleem van Sturm-Liouville, Schrijf (5.25) 
als volgt: 


Lw + Apw= 0; Les + Q(x). (5.24) 


Definitie (5.4): 

Een differentiaaloperator Ì heet zelf-geadjungeerd indien geldt: 
uiv-viu= Et f(u,v‚u',v') 

waarin f een geschikt gekozen funktie is, 


Operator L (5.24) is zelf=-geadjungeerd immers: 
uLv-vLu= En [uv' — vu!) . (5.25) 


Als geldt: u(O)=u(1)=v(O)=v(1)=0, dan wordt (5.25) na integratie: 
kl 
| (uL v-vLu)dx= 0, (5,26) 


Voor het probleem van Sturm-Liouville (5,23) gelden nu de volgende stellingen: 


Stelling (5,2):(zonder bewijs) 


Bij iedere eigenwaarde An behoort één en slechts één eigenfunktie Dn 


Stelling (Ba): 


De eigenwaarden A en bijbehorende eigenfunkties Ge zijn alle reëel. 


Bewijs: 


De operator L is reëel, p(x) is een reële funktie, Als de complexe A en P, 
aan het probleem voldoen, dan voldoen ook de toegevocgd complexe Le en 9 
(ka na). Dan geldt: 

L(P) = =APP, 

L(9) 
We berekenen SG) (Dax Ag) fo lolfax. 
Met behulp van (5.26) concluderen we dat het linkerlid nul is, Dus 


4 
(An … 5, | PP Pp dâx= 0. 


-1,09,° 
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EN 


Daar Pov =elo, | definiet positief is op het interval, geldt: 
AA 

Dus is A, reëel. 

stel bij een reële eigenwaarde An hoort een complexe eigenfunktie De 

Dan voldoen An en de toegevoegd complexe ®, Eveneens aan het problecen, 


Volgens stelling (5.2) geldt dan 00 Dus op, is reëel, 


Stelling (5.4): 


De eigenfunkties ’, van het probleem voldoen aan de orthogonalitcitsrelatic 





4 
dE o(x) PPS = 0 indien n + m 


f 
p(x) pApdx FO indien n = m. 


cd 


als Pen P. eigenfunkties zijn van het probleem, geldt: 


Lia) = =A, 6 9 
Lo) 


We berekenen: 


ai p Pin 


1 
f Lonb(o,) — oblo) dxe= (A, - A) P 0, 9% dt. | 
od à 


Met behulp van (5,26) concluderen we dat het linkerlid nul is. Dus 
f 
On - A) d P 9, PM dx = 0, 


Als m *n geldt de orthogonaliteitsrelatie 


4 
on Po % dx=0, (BeB7) 


We zullen in het volgende veronderstellen dat de cigenfunkties zodanig genor= 
meerd zijn, dat geldt: 
\| 
Ì Pe, ded, | (5.28) 
o 


De eigenfunkties vormen een orthonormaal stelsel van funktics, 


Opmerking: 
Zowel in stelling (5.3) als (5.4) is gebruik gemaakt van de ze 





heid van de operator L, 
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stelling 1Baide 


Voor de eigenwaarden en (n=1,2,...) van het Sturm-Liouville probleem bestaat 


een retel getal a zodanig dat A >a (nld pende 


Bewijs: 
Vermenigvuldig (5,23) met w en integreer over het interval. Dit levert: 


1 2 d | 
| wW 2E ax + J Md + nf p a = 0, 
dx 


Partiële integratie van de eerste term levert met behulp van de randvoorwaarden: 


1 1 rl 
| p nik JE Qw'ax = d (2) dax. (5.29) 


Nu geldt ej p wd >» O0 en Er À (Ee dx > O. Hieruit volgt: 


a) Als Q definiet positief, dan geldt A> O0, 

b) Q(x) kan negatieve waarden aannemen op het interval, maar |e(x)| en 
la(x)| zijn begrensd op [0,1]. 
Daar p>0 kunnen we altijd een getal ad vinden zodat 


„ap + Q < O0 voor alle 0 Sx S 1, 


We stellen A=Ì-a, dan wordt (5,29): 


| Á | 
| pan | ey? a | lapeëtle® a, (5.30) 


Dan geldt ADO en dus AD-a, 
c) Het geval dat lo | en of |a| niet begrensd zijn op [o,1] behandelen we niet. 


Enkele algemene resultaten uit de theorie van Sturm=Liouville, 


We schrijven vergelijking (5.24) in de vorm 


Llw) = F (Be A1) 
me t | 

F = =Àpw (5,32) 
met randvoorwaarden: 

w(O) = w(1) = 0, | (5.53) 


Op vergelijking (5,31) met randvoorwaarden (5.33) kunnen we formeel formule 
(5,12) toepassen. Voorwaarde hiervoor was, dat de homogene vergelijking 
L(w)=0 geen oplossing had die aan de randvoorwaarden voldeed. Dit komt overeen 


met de eis dat A=0 geen eigenwaarde is van het SturmeLiouville probleem, 


Hulpstellin „6): 


Deze eis is geen beperking van algemeenheid, 


60, 


Bewijs: 


Stel ÀA=0 is een eigenwaarde van het probleem. De konstante u Komt niet voor 
onder de eigenwaarden van het probleem, Tel bij het linker- en rechterlid van 


vel. (5.23) op bew en stel: 
L= L + up (5.34) 


>et 


= el + À 
dan worden vgl. (5.31) en (5,32) gekombineerd 

L(w) = =Aow (5e) 
met de randvoorwaarden (5,33) terwijl A=0 geen eigenwaarde van het probleem is, 


Wij nemen vervolgens aan dat na transformatie het probleem geen eigenwaarde 


A\=0 toelaat, Wij schrijven de oplossing van (5,31) in de vorm: 


1 
wlx) = | P(E)G(x, E)dE (5.36) 
of met (5,32) 


1 
w(x) = =A d w(E)o(E)G(x, E)dE, (5.5) 


We merken op dat w(x) een oplossing is van de integraalvergelijking (5.37). 
Voor dit type vergelijkingen bestaat een uitvoerige theorie die ons in staat 
stelt een aantal fundamentele eigenschappen van het probleem Sturm-Liouville 
te bewijgene. 


Zonder bewijs vermelden we: 


De existentiestelling . 

Er bestaat een oneindige rij van reële eigenwaarden de met de eigenschap 

A1 Àn ‚ n=1,2, Sg eee 

Als Q(x)>0, zijn alle eigenwaarden positief; is Q(x)<O dan bestaan ven eindig 
aantal negatieve cigenwaarden, 


A, neemt onbegrensd toe als noo, 


De ontwikkelingsstelling (5,8): 

Zijn P de benne eigenfunkties van het Sturm=Liouville probleem, dan kan 
elke funktie vx) ec? [0,1] die voldoet aan de randvoorwaarden (5.33) ontwikkeld 
worden in een op [O,1] uniform en absoluut konvergente rveks van de vorm: 


oo 
Y(x) = DN À, 9 (x) 
n=0 


| 
A= | PP, v(x)ax, 


met 
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Toepassing: ontwikkeling van de funktie van Green. 


Beschouw de inhomogene vergelijking met randvoorwaarden 


afs 
S5 + [a(x) + Aolx)]w = F(x) (5,38) 
dx 
w(O) = w(1) = 0 (5.39) 


met gegeven funkties Q(x), e(x), r(x)ec [0,1]; A is konstant, p>0 op (0,1). 
Stel A is geen eigenwaarde van het homogene probleem (probleem (5.38) en. (5.39) 
met F(x)=0). 

De oplossing w(x) van het inhomogene probleem bepalen we nu met behulp van de 
ontwikkelingsstelling als reeksontwikkeling, waarbij de termen op een konstante 
na gevormd worden door de eigenfunkties van het homogene problcem, 

Als ’, een eigenfunktie is en A de bijbehorende eigenwaarde van het homogene 


probleem, dan vindt men met eenzelfde redenering als bij eigenschappen 2 en 5: 


4 1 
a-) | en va | °F äx. (5,40) 


We voeren in 


oo 
w(x) = E b, Dy 0). (5.41) 
k= 
Vel. (5,41) met PP, vermenigvuldigen en integreren levert voor n=1,2, 54 «ee 
1 
pel pF dx 


bo = (5.42) 


n EN : 
n 
In plaats van de oplossing te ontwikkelen naar eigenfunkties kunnen we met 


behulp van de theorie van de Greense funktie de oplossing als volgt schrijven 
1 É 
wa) si C(x, E)F(B)äE. (5.43) 


Uit (5,41) — (5,43) volgt: 
co glx), (5) 


G(x,E) = DE an d 
= n 


í 


Met andere woorden, de funktie van Green voor het inhomogene probleem hebben 


n 


we gekonstrueerd in de vorm van een reeks, waarbij de termen bepaald worden 


door de eigenfunkties van het homogene probleem, 


We veronderstellen dat A in (5,38) geen eigenwaarde was van het homogene 
probleem, Stel daarentegen nu dat A=À (p een natuurlijk getal). We schrijven 


de gekonstrueerde oplossing (5.41), (5,42) in de vorm: 


1 Ge), (5) a) ff 
dl EN 2 IR Dat + Be | sr Hat. 
n n Pp 
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De oplossing bestaat dus uitsluitend in deze vorm als 


4 
dl (OF (Hat = 0 


In dit geval is de oplossing niet volledig, Wij schrijven de oplossing dan 


in de vorm: 


1 9), (6) 
w(x) = | sed ne EE + an % 


AA 
n 


waarin a een konstante is, 


0. 
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Hoofdstuk 6: Vraagstukken, 
Hoofdstuk 1: 
Welke van de volgende differentiaalvergelijkingen zijn lineair? Welke van de 


lineaire differentiaalvergelijkingen zijn homogeen? 
2 


a) 2 +r=e 
de 
dy e_ 1 
Oe ae 


Xx 


2 
c) AX + (1 + cosx)y = sinx 


dx 
2 
ä) x° sE + x 2e + (x= 1)y = 0 
dx 
2 
d 2, d 
jep ler) Ee 
dx 


ae 
£) se + siny = 0, 
dx 


Bepaal ook de orde van bovenstaande differentiaalvergelijkingen,. 
Schrijf 1 a, c,‚ d, e, f als een stelsel van eerste orde differentiaalvergelij=- 


kingen. 


Bepaal de algemene oplossing van: 


2 
SE 2 dy en zi = SinXe 
dx 


welke van de volgende funkties 


Ya = xsinxs Ya = Reede : Yz = Yx COSX 
1 2 5 


Vx 


is een oplossing vans: zy + XY! + fs = Dr = 0, 


Gegeven zijn de differentiaalvergelijkingen Epi = 2Yy en x = In(1+sin se In 
welke punten is Vy en In(1+sin zij niet Lipschitz kontinu? 


Ga na wat dit betekent voor ed = 2Vy y(o) = 0, 


Hoofdstuk 2: 


Bepaal het richtingsveld van de differentiaalvergelijkingen: 


a) 2 en kg 


ï 
ee 
Be 


d 2 2 
b) De = 


8. 
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11, 


14. 


15 


Bepaal de algemene 


2 2 
(3x°-y)y' = 2xy, 


64, 
oplossing van de vergelijkingen 
U (1-y). 


dx 


y' = 5 (1 + Iny =— Inx). 


oplossing van ay = 


Bepaal de algemene Er 
vlak, 


XY 


X+y en geef een schets van de integraal 


krommen in het 


Bepaal de algemene oplossing van Ee + y = SinXe 


LE a dmt beginvoorwaarde y(O)=1. 


Bepaal de oplossing van dx or 


Bepaal de algemene oplossing van: 


2d è 
Xx De = Xy + y = 0 


x 


dx 0, 


2 
-4y - Xx V\y = 
Toon aan dat de volgende differentiaalvergelijkingen exact zijn en bepaal de 
algemene oplossingen: 


2 dy 


2 
Ba & 3x Y + 2xy = 0 


a x° + 2y) 


(x-cosxy + e)) Een + xy cos xy + sin xy = O0, 


Bepaal alle differentiaalvergelijkingen, van het type 

d 

Ee = f(x,y) 
die een integrerende factor M=X(x).Y(y) hebben. (X(x) is een niet-konstante 
funktie die alleen van x afhangt en Y(y) is een niet-konstante funktie die 


alleen van y afhangt.) 


Bepaal de integrerende faktor van: 


a + glx)y = h(x). 


Bepaal vervolgens de algemene oplossing en vergelijk deze met formule (2,22), 


diktaat blz. 6, 


PES 
De 


17, 


18 


20. 
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Gegeven: 


dr 2 a 
a ay + bx a,b £ O, 


2) Onderzoek de gevallen a=0 en a==2, (Stel voor dit laatste geval y-l) 
b) a tE -j, 
Beschouw de transformatie: 


a+ Î 
xX 





sk 
a+1 y 
Bepaal de differentiaalvergelijking voor Y=Y(X). 


c) Beschouw de transformatie: 


0 pn 1 Eed 

Kep Pet En). 
Bepaal de differentiaalvergelijking voor n=n(E). 

d) Leidt uit b) en c) af voor welke waarden van a de oorspronkelijke differentie 


aalvergelijking elementair oplosbaar is, 


Bewijs dat de orde de differentiaalvergelijking Blaren Peo die homo= 
geen is van de graad k met één verlaagd kan worden door te stellen 
y(x)=exp(lu(x)). De differentiaalvergelijking Player” sonen 22=0 heet homogeen 
van de graad k als geldt: 


(n) 


k 
Bl der hart aen re 3 = A Flusyor oneens? 


. Ì 
\ 


Bepaal de algemene oplossing van: 


yy - (y=xy')Ê = 0 (maak gebruik van vraag 17). 


Gegeven is dat ay = f(x,y) een integrerende faktor M=M(x) heeft. 


Bewijs dat de differentiaalvergelijking lineair is. 


Bepaal de integrerende faktor van: 


2xy + xy + E y + (x2 ey) En = 0, 


Is dit antwoord in overeenstemming met 19? 


Bepaal de algemene oplossing van: 


dy _ 
dn y tg Xx + COSX. 


Hoofdstuk 5: 


2) Schrijf ja, c, d, e, f‚ als een eerste-orde vector differentiaalvergelijking. 


_ 


b) Bepaal voor la en 1e y(x), g(x) en Alx) in: 1 = Alx)y+ g(x). 
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23. Gegevem de differentiaalvergelijking: 


dy zr L > en 
Te ° f(y‚x) en beginvoorwaarde y(x) Jo 
y(x) en f(y‚x) zijn n-dimensionale vectorfunkties. 

In de (n+1) dimensionale Euclidische ruimte is het gebied G gedefinieerd doors 


[y-y, | = b, [xx | sa (a en b zijn konstanten), 


Verder geldt: 
EG) sk) 1 + {71} 


EGsx) - Er) sk) 17,-9| 


HA 


Ys VEE Jo € Ge 
k(x) is een kontinue funktie in x. 
f(y‚x) is kontinu op G in x en y. 


bewijs met behulp van de methode van de suksessieve approximaties de existentie 
van een oplossing y(x) van de differentiaalvergelijking, waarvoor geldt 


vC )=T Bewijs dat voor deze oplossing geldt: 
vG) - Yol SI + Iv, |) (exo( je k(x')dx') = 1) 
x . 
o 


voor [x-x | = a 
o | 


a is een konstante waarvoor geldt «sa, 


24, Bepaal de oplossing van: 


dy _ = 
rek A: y(o) = 0 
met behulp van de methode van de suksessieve approximaties, 


Vergelijk het resultaat met het antwoord van vraagstuk 9, 


25. Gegeven de vector-differentiaalvergelijking: 


_ ’ 


WY Ay + F(x) 


dx 
waarin: 
y an (74192) 
0 1 
ais G 0) 
F(x) = (0,x) 
y(0) zl Yo ad (Ce, sen) 


e“ en e* zijn twee lineair onafhankelijke oplossingen van het homogene stelsel 


differentiaalvergelijkingen, 


zip 


27. 


23, 


350, 
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a) Bereken Ú(x) in Yi = EPA waarin ke oplossing is van het homogene stelsel 
en P(x) een matrix waarvoor geldt: f(O)=I (I=eenheidsmatrix). 
b) Bereken de algemene oplossing y van het inhomogene stelsel. 


c) Bereken de oplossing van het inhomogene stelsel als Yo=(240). 


Gegeven de vector differentiaalvergelijking: 
dy rn 
en f(y‚x). 

yen f 

Onder welke voorwaarden kan men uit dit stelsel de volgende differentiaal- 


zijn n-dimensionale vectoren met komponenten Y; resp. fj. 


vergelijking afleiden: 





n-=1 
n e5 NEER, ) 
dx 


Gegeven is het stelsel: 


dy dy dy 
ds KA nn rs 
de Te a NE TS Mt Ape 


Bepaal de algemene oplossing, 


Gegeven is de vector differentiaalvergeiijking: 


e ei Á -@ 
d _— — 

SL Ay met Y = (y4Y5sYz) en A= ( O =1 4) 
dx tha Aa 4 O0 4 


Bepaal de algemene oplossing, 


Bepaal de algemene oplossing van SL = Ay waarin: 


== ú % 
A= Kh 6) en ye (y45Y2)e 


Hoofdstuk 4: 


Gegeven is de derde orde lineaire inhomogene differentiaalvergelijkings 
y"' pen 2y“ pen y! + 2y = 4Xe 


Bepaal de algemene oplossing, 
schrijf de vergelijking in veetorvorm, bereken de fundamentele matrix. 


Bepral de oplossing waarvoor geldt: 


y(0) = 0; y'(0) = 2; y"(0) = -1. 


Sha 


Ss 


54. 


55e 


68, 


Bepaal de algemene oplossing van: 


Ee 2 } 
5 + Wo x= Â sin wt. 
dt 


Eoe luidt de oplossing als wew ? 
Ca de voorwrarden van de existentie= en eenduidigheidsstelling na voor: 
ex yr — 2xy! + n(n+1)y = O0 __n > O0, geheel. 


Bepaal voor n=1 de algemene oplossing. 


Benader de oplossing van: 


dx 


met behulp van; 


ä 2 1 
Lex-y y(o) = 5 


le 5 (xx) dy(x) 
Ws Bme 
ie n=1 a! dx 


ij 


a) y(x) 


b) y(x) 


u 
® 
+ 
u a 
@® 
al 
Ld 


Herleid een der oplossingen tot een exponentiële funktie, 


Gegeven de differentiaalvergelijking: 





Ly in f(x), 
me t tiÂ 
d d. d 
L= a == + a + eee A, TT + A 
n n n= n-Î 1 dx o 
dx dx, 


waarin anr*"1äg konstanten zijn. 


Ly=0 heeft eigenwaarden As i=Îgeees Ne 
Wij onderscheiden voor de bepaling van een particuliere oplossing Yp van de 
inhomogene vergelijking de volgende gevallen: 


n) f(x) =a el; a=konstante; A, us dief an on ole 
Stel adi ae 


Voorbeeld: 


gt „yy! +y! y= e°*, 


57e 
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b) f(x)=a e*; n= A, terwijl A, een multipliciteit s heeft, 


Stel Jy =P, me * x° waarin Fitt) een veelterm van de graad s-1 is, 
Voorbeeld: 
y' + 27! +y=e%, 
c) f(x) =P (2); Stel vpn) P(X) en ax) zijn veeltermen van de graad me 


Voorbeeld: 


„4 5 4 


+ 2y! + y= XxX + 2x = sr + 1e 


d) f(x)eP, (jet Agtus iel,eesjne 
8 ux 
Stel Toa de : 


Voorbeeld: 
y' + 2y! + y= X e”“, 


e) £(x)=P (oe %; hel, terwijl A een multipliciteit s heeft. 
n _ Ss ux 
Stel Jp An) x ee”, 


Voorbeeld: 
_X 


a A Et 


Bepaal met behulp van machtreeksontwikkelingen: 


a) de algemene oplossing van: y'+(sinx)y=e" 3 


b) een particuliere oplossing van y"+2xy'+4y=e". 


Voor welke waarde van ad heeft 


y!" = 2xy!' + ay = 0 


een oplossing die een veelterm is in x? 


Gegeven de differentiaalvergelijkingen: 
2 


a) a + sinp = 0 
dt 
a Î 
b) an ik ad Âe 0 
dt 
de 
ec) Sa + p= 0, 
dt 


Voor alle vergelijkingen gelden de beginvoorwaarden: 


(0) =o LO) = 


Bepaal de oplossing van a) door de eenen agnes achtereenvolgens 


te differentiëren en de verkregen waarde van a’ (o) te substitueren in de Taylor= 


reeks van p(t). at” 
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41. 


42, 


45. 
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Schrijf b) als een stelsel en bepaal met behulp van de methode van de sukses- 


sieve approximaties de oplossing, 


Substitueren in c) pet, bepaal A; ontwikkel de verkregen oplossing in een 


Taylorreeks, 
Vergelijk de coëfficienten van de reeksontwikkelingen van de oplossingen van 


de drie differentiaalvergelijkingen,. 


Bepaal de algemene oplossing van: 


2 
5 ar - 2y = pd + RL voor x > 0, 
dx ze 


Bepaal de index-vergelijking en de recurrente betrekkingen voor de coëfficienten 


van de machtreeksen, die oplossingen zijn van: 


a) 2xy' + y' + xy = 0 (x > 0) 
b) xy" +y' -y = 0 (x > 0), 


Bepaal door middel van machtreeksontwikkelingen twee lineair onafhankelijke 


oplossingen van: 
Ey + XY! + (x°-1)y = 0 * >.0, 


Ga na waarom ln x niet in de algemene oplossing voorkomt. 


Bepaal van de differentiaalvergelijking van Laguerre: 
xy! + (lex)y' + Ay = O Ca > Oy 


(A is een positieve konstante), 
De index-vergelijking, de recurrente betrekking en een oplossing. Toon aan dat 
als À een positief geheel getal is, deze oplossing zich reduceert tot een veel- 


term. 


Hoofdstuk 5: 


Gegeven is de differentiaalvergelijking 


2 
y' + 5 y= f(x) 
met f(x)€C'[0,1] 


en randvoorwaarden y(O)=y(1)=0. 
Bepaal de funktie van Green voor dit probleem en geef de oplossinge 


44. Gegeven is de differentiaalvergelijking: 


xy + xy! -— e y= Yx O<x5sie (1) 


45. 


46, 


41. 
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a) Breng de homogene vergelijking-m.b.v. de transformatie van Liouville in de 
standaardvorm, 
Bepaal vervolgens twee lineair onafhankelijke oplossingen u, en u, van de 
homogene differentiaalvergelijking in de standaardvorm zodanig dat: 
u, (1) = 0; lim u(x) = 0. 
xJo 
Bepaal de funktie van Green. 
b) Bereken y(&), als y(x) voldoet aan (1) en aan de randvoorwaarden y(1)=0; 
lim y(x)=0. 


XO 


Gegeven het randwaardeprobleem: 


WA 
ed 
A 
® 


y' + 2y! + y= f(x); 1 


y(i)=yle)=0 F(x) EC'[1,e). ne 


Bepaal de funktie van Green voor dit probleem; bereken y(2) els f(x)= 5e 
x 





Gegeven: 
è 


ar + In = sinkx 
dx 


y(0) = y(r) = 0 
A en k zijn positieve konstanten. 
a) Bepaal de oplossing van dit probleem met behulp van de funktie van Green, 
b) Geef de ontwikkeling van de funktie van Green in de eigenfunkties van het 
homogene probleem, 
c) Zij A een eigenwaarde van het homogene probleem, Voor welke waarde(n) van 


k bestaat dan de oplossing? 


Gegeven: 
2 
Sa ak er 
2 ä x dx 2 y f(x) 
dx xX 


waarin f(x)€C'[1,4] 
en randvoorwaarden y(1)=y(4)=0, 
Bepaal de oplossing van dit randwaarde probleem met behulp van de funktie van 


Green. 


72. 


48, Gegeven: 
2 
Ë ai — 2y = f(x) 
dx 


f(x) € c'[o,1)s (0) = 0; YI) = a, 
a) Bepaal als a=0 de oplossing van het randwaardeprobleem met behulp van de 


funktie van Green, 


b) Bepaal de oplossing voor het geval dat a=1, 
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Het totale werk is 20 pt. waard. 


Opg.1 = 5 pt. 
Beschouw voor t > 0 het stelsel diff. vgl. + beginvoorwaarden 


(ed. 


EE Ze exp(-y) s x(0) = 0 
(*) 
ay = exp(-x) 3 y(O) = 0 


a. We definiëren de baankromme behorend bij (*) als de deel- 
verzameling van het [ x,y l-vlak van de punten [ x(t),y(t) | met 
t > 0 en x(t), y(t) de oplossingen van (*). 


Zij deze baankromme gegeven, als y = f(x) voor x= 0. 


Laat zien, dat de helling van de baankromme dy voldoet aan 


dy … dy/at 
dx dx/dt" 


b. Geef een diff. vgl. + beginvoorwaarde voor de baankromme en 
los deze op. 
Geef een schets van de baan. 


Bestaat er een asymptoot voor t >» ? Zo ja, welke ? 


c. Bepaal tenslotte x(t). 


Opg.2 = 5 pt. 


Gegeven zij het onderstaande inhomogene stelsel diff. vgl. 


el Oe een 


a. Bepaal 2 lineair onafhankelijke oplossingen van het homogene 


be han 
ct 


probleem. 


‚0. 


b. Geef vervolgens de algemene oplossing van het inhomogene 


probleem. 


opg. 3 = 6 pt. 
Beschouw de volgende integraalvgl. voor t > 0 


s(t) 
(*1) y(t) = Yo + S f(Tsy(T))drt met yo ER. 
0 


Er geldt 
(1) sE cit o „al s s(0) = O0 en voor 0 <t <a: 0Ss(t) St 


(ii) f E C(G) en f(t‚y) is Lipschitz-continu in y met Lipschitz- 
constante A op het gebied G met G = [O,al Xx [yobs Yo*b] 
met a,b > 0 


Beschouw naast (#1) de diff. vgl. + beginvoorwaarde 
(#2) y(t) = sT(t). fltyyls(t))) 3 y(0) = Yp: 


a. Laat zien,dat geldt 
ye cht o ‚al en y is oplossing van (*2) © 


yEecCf[O0,al en y is oplossing van (*1) 


Definieer in dit geval een Picard-Lindelöf iteratieproces. 


4 


ec. Laat zien, dat dit proces uniform convergeert in C[ Gol 
met geschikte a > 0 naar een oplossing van (*1) 
Hint : Vind een schatting voor het verschil tussen twee op- 


eenvolgende iteranden en gebruik de idenditeit 


n 
z Yo + 2 (y‚-y;_4) om aan te tonen, dat y_ een 
0 iz È 7E -, n 


Cauchy rij is. 


In 


Opg.b = 4 pt. 


Als een speciaal geval van opg. 3, (+2) beschouwen we 


z {td = zkt) # 20) = d 





| 
\ 


oo 
a. Geef de oplossing in de vorm van een machtreeks z = & an t 
waarbij de coëfficienten a, expliciet zijn bepaald. 
b. Het is interessant het groeigedrag voor t > @ van z te 


vergelijken met exponentieel groeigedrag voor t > ©. 
Toon aan, dat voor. elke à > 0 geldt 


z(t) = olexplAt)) voor t > o, 


" n 
Je © f Gek aa ie VEE 
an gl bo. eelan den Exite sap. ORN en 
Sind, dre Îv [x, d, enn la ba den Gdret 

S jd, Ja) fig „Je pt olurch j elen 

) 

Plt P ct olen Kon dina ben 8, Nn -- W 

SON 4 minder lems Cing Jk legal kun VE, nd 
dizst barn had Lerdo, Sette. bin an der 

el a N forgoeht werde, 

Ut: pi eran he alleich van 

ones id ole Urd LS tn gen 
Kan he 


A Cd 2-l, 


Tentamen Grondslagen Differentiaalvergelijkingen 28 mei 1974, 
13.30 = 16.30 uur. 


Beschouw de vergelijking van het type Euler 


tt tt 


3 2 ' R 
Xx y + 2X y + AXy =— Ay = 0 
voor Xx > 0 5 A is een willekeurige konstante. 
Bepaal drie onafhankelijke oplossingen van de vergelijking. 


Onderscheid hierbij verschillende waarden van À en bewijs de 


onafhankelijkheid van de verkregen oplossingen. 


Gegeven is de vergelijking 
ge + Ay = 0 
De konstante À wordt bepaald door nader te speciferen voorwaarden. 
a. Beschouw de vergelijking op het interval [0,1] met randvoor- 
waarden y(0) = y(1) = 0. 
Geef de eigenfuncties en eigenwaarden van het probleem. 


b. Beschouw het interval [1,3] met randvoorwaarden y(1) = y(3) = 0. 


Geef opnieuw de eigenfuncties en eigenwaarden. 


Beschouw de vergelijking 
nk n uy & 2ay = 0 
a is een konstante. 
a. Bepaal twee lineair onafhankelijke oplossingen door reeks- 
ontwikkelingen in een omgeving van x = 0. 
We noemen deze h‚ (xx) en h(x). 
Bepaal de konvergentiestraal van de reeksen. 
N.B. In de volgende onderdelen spelen h‚ 60) en hx) opnieuw 


“een rol; het is echter niet de bedoeling dat hierbij de 


verkregen reeksontwikkelingen worden gesubstitueerd. 


b. Beschouw vervolgens het randwaardeprobleem 


tt 


yy - Zig + 2ay = fx) 
met randvoorwaarden y(-1) = y(1) = 0 5 
xEl-1,1l en f(x) € H-1,1l; f(-1) = £(1) = 0 
De konstante a is zodanig gekozen dat er geen niet-triviale 
oplossingen van het homogene randwaardeprobleem bestaan. 
Breng door middel van de Liouville-transformatie het probleem 


in de gedaante 


tt 


nt RGoOn = £*(x), nl-1) = nl1) = 0 (+) 


Bereken de wronskiaan w behorende bij het oorspronkelijke 

en het getransformeerde homogene probleem. 

Indien n, en nj) lineair onafhankelijke oplossingen zijn 
van het getransformeerde homogene probleem, luidt de oplossing 


van het inhomogene probleem (zonder toepassing van randvoor- 
waarden) 

Ea € 
RAS ed En, Ge), (ED = ngen, (EN £ (BAE + An, Go) + Ayn, Go) 


Bepaal de Greense functie behorende bij het randwaardeprobleem 
en geef een uitdrukking voor de oplossing y(x) van het oor- 


spronkelijke randwaardeprobleem. 








